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A model of the passive vector quantity advected by the Gaussian velocity field with the covariance 
oc S{t — t')\x — x'l^ is studied; the efi'ects of pressure and large-scale anisotropy are discussed. The 
inertial-range behavior of the pair correlation function is described by an infinite family of scaling 
exponents, which satisfy exact transcendental equations derived explicitly in d dimensions by means 
of the functional techniques. The exponents are organized in a hierarchical order according to their 
degree of anisotropy, with the spectrum unbounded from above and the leading (minimal) exponent 
coming from the isotropic sector. This picture extends to higher-order correlation functions. Like 
' in the scalar model, the second-order structure function appears nonanomalous and is described 

' by the simple dimensional exponent: S2 oc r^~'^. For the higher-order structure functions, S2n oc 

e)+ „ ^ ^j^g anomalous scaling behavior is established as a consequence of the existence in the 
5— ( ' corresponding operator product expansions of "dangerous" composite operators, whose negative 

critical dimensions determine the anomalous exponents A„ < 0. A close formal resemblance of the 
model with the stirred Navier-Stokes equation reveals itself in the mixing of relevant operators and 
is the main motivation of the paper. Using the renormalization group, the anomalous exponents are 
calculated in the 0{e) approximation, in large d dimensions, for the even structure functions up to 
the twelfth order. 
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I. INTRODUCTION 



The investigation of intermittency and anomalous scaling in fully developed turbulence remains essentially an open 
^ ' theoretical problem. Much effort has been invested recently into the understanding of the inertial-range behavior of 
04 the passive scalar. Both the real experiments and numerical simulations suggest that the breakdown of the classical 
' Kolmogorov-Obukhov theory is even more strongly pronounced for a passively advected scalar field than for the 
, turbulent velocity itself. On the other hand, the problem of passive advection appears easier tractable theoretically; 
see Ref. [|| and references therein. 

The most progress has been achieved for the so-called rapid-change model of the passive scalar advection by a 
■ self-similar white-in-time velocity field The model is interesting because of the insight it offers into the origin 
of intermittency and anomalous scaling in turbulence: for the first time, anomalous exponents have been calculated 
on the basis of a microscopic model and within controlled approximations 0-0]. Within the "zero-mode approach" 
to the rapid-change model, proposed in Refs. nontrivial anomalous exponents are related to the zero modes 

(homogeneous solutions) of the closed exact differential equations satisfied by the equal-time correlations. In this 
^ ' sense, the model appears exactly solvable. A recent review and more references can be found in Ref. @. 
k>( I In Ref. 1^ and subsequent papers the field theoretic renormalization group (RG) and the operator product 

j_] ■ expansion (OPE) were applied to the model The feature specific to the theory of turbulence is the existence 

d ' in the corresponding field theoretical models of the composite operators with negative scaling (critical) dimensions. 
Such operators, termed "dangerous" in |p[-p^, give rise to anomalous scaling, i.e., the singular dependence on the 
infrared (IR) scale with certain nonlinear anomalous exponents. 

The OPE and the concept of dangerous operators for the Navier-Stokes (NS) turbulence were introduced in Ref. [ p5[ ; 
detailed review and bibliography can be found in p^ , |r7[ . The relationship between the anomalous exponents and 
dimensions of composite operators was anticipated in Ref. for the stochastic hydrodynamics and in |0,^,|l9| for 
the Kraichnan model within certain phenomenological formulation of the OPE, the so-called "additive fusion rules," 
typical to the models with multifractal behavior |po| . A similar picture naturally arises within the context of the 
Burgers turbulence and growth phenomena [ pl|j2^ ]. 

Important advantages of the RG approach are its universality and calculational efficiency: a regular systematic 
perturbation expansion for the anomalous exponents was constructed, similar to the well-known e-expansion in the 
theory of phase transitions, and the exponents were calculated in the second [p|~|rT|] and third orders of that 
expansion. For passively advected vector fields, any calculation of the exponents for higher-order correlations calls for 
the RG techniques already in the 0(e) approximation Furthermore, the RG approach is not related to the 
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aforementioned solvability of the rapid-change model and can also be applied to the case of finite correlation time or 
non-Gaussian advecting field jlj] . 

Recent research on the Kraichnan model and its descendants has mostly been concentrated on the passive scalar 
advection. The large-scale transport of vector quantities exhibits more interesting behavior; see monograph and 
references therein. In this paper, we study the anomalous scaling and effects of anisotropy and pressure, for the 
passive vector field advected by the rapid-change velocity field. The model has already been introduced and discussed 
independently in Rcfs. and 

Before explaining our motivations, which follow the same lines as those of Refs. |2^,^, we shall discuss the definition 
of the model in detail. 

We shall confine ourselves with the case of transverse (solenoidal) passive 9{x) = {9i{t,x)} and advecting v(x) = 
{vi{t,x.)} vector fields and the advection-diffusion equation of the form 

Vt0^ + ^,r = :^oM, + f„ Vt=dt + iv,d,), (l.I) 

where V{x) is the pressure, z/q is the diffusivity, A is the Laplace operator and fi{x) is a transverse Gaussian stirring 
force with zero mean and covariance 



{Mx)f,{x')) ^6{t-t')aj{r/L), r = x-x'. (1.2) 

The parameter L is an integral scale related to the stirring, and Cy is a dimensionless function finite as L — > oo. 
Its precise form is not essential; for generality, it is not assumed to be isotropic. Therefore, the force maintains the 
steady state and is also a source of the large-scale anisotropy in the system. 

The velocity v{x) obeys a Gaussian distribution with zero mean and covariance 

{v,{x)v,{x')) = Do5{t-t') JvpP,j{p)p-''-' exp [i(pr)] . (1.3) 

Here and below p is the momentum, p = |p|. Dp — dp/{2TT)'^, Pij{p) = Sij — PiPj/p^ is the transverse projector, 
Do > is an amplitude factor and d is the dimensionality of the x space. The exponent < £ < 2 plays in the RG 
approach the same role as the parameter e = A — d does in the RG theory of critical behavior. The IR regularization 
is provided by the cut-off in integral ( |1.3| ) from below at p = m, where 1/m is another integral scale; the precise form 
of the cut-off is not essential. In what follows, we shall not distinguish the two IR scales, setting m ^ 1/L. The 
relations 



Do/ 1^0 = .90 = 



(1.4) 



define the coupling constant go (i.e., the formal expansion parameter in the ordinary perturbation theory) and the 
characteristic ultraviolet (UV) momentum scale A. 

Due to the transversality conditions, diOi = diVi = 0, the pressure can be expressed as the solution of the Poisson 
equation. 



AT' = -d^vjdj9^. 

The issue of interest is, in particular, the behavior of the equal-time structure functions 



(1.5) 



(1.6) 



in the inertial range, specified by the inequalities 1/A ^ r <C i ~ 1/m. Here 9r = 9iri/r is the component of the 
passive field along the direction r = x — x', an analog of the streamwise component of the turbulent velocity field in 
real experiments. 

The general symmetry of the vector problem permits one to add on the left-hand side of the advection-diffusion 
equation the "stretching term" of the form {9jdj)vi. This general vector model is studied in Rcf. ||2^, and its special 
magnetic case, where the pressure term disappears, was studied earlier in a number of papers in detail; see Refs. 

From the physics viewpoints, the model (1.1 )-( pT3| ) can be considered as the linearized NS equation with the 
prescribed statistics of the background field v and the additional convention that the field v is "soft" and the 
pertur batio n 6 is "hard," that is, 89 ^ dv. Our motivation, however, is different: a close formal resemblance of the 
model (1.1)-(1.3) with the stirred NS turbulence. 

Although the model (1.1) is formally a special case of the general "y^-model" ||2^, it appears in a sense exceptional 
and requires special investigation. In this case, the stretching term is absent and the analog of the kinetic energy 
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9'^{x) is conserved, as for the passive scalar and stochastic NS equations. An important consequence is the existence 
of a constant flux solution, characteristic of the real NS turbulence 0] . The second-order structure function appears 
nonanomalous and is described by the simple "dimensional" exponen t, Sp , oc r^"*^. 

Furthermore, since only derivatives of the field 9 enter into Eq. (1.1), the latter possesses additional symmetry 
6 —I- 9 + const. The leading anomalous exponents in the magnetic and general cases are determined by the composite 
operators built of the field 9 without derivatives in our case they become trivial and the leading nontrivial 

exponents are related to the composite operators built solely of the gradients of 9. 

A similar distinction exists also between the density and tracer scalar fields advected by a compressible velocity, 
as discussed in Ref. Q in detail. But in contrast with the scalar case, where the leading exponent for a given-order 
structure function is determined by an individual composite operator even in the presence of the 9 9 + const 
invariance ||^ in the vector model (1.1) the inertial-range behavior of any structure function is determined by a 
family of composite operators with the same symmetry and dimension, and in order to find the corresponding set of 
exponents and to identify the leading contribution, one has to consider the renormalization of the whole family, which 
implies the mixing of individual operators. 

In the scalar case, the anomalous exponents for all structure functions are given by a single expression which 
includes n, the order of a function, as a parameter |^,^j8|. This remains true for the vector models with the stretching 
term p3| , ^ . In the special vector model (IT), the number and the form of the operators entering into the relevant 
family depend essentially on n, and different structure functions should be studied separately. As a result, no general 
expression valid for all n exists in the model, in contrast with the scalar |^,|^, magnetic and general vector ]2^ ] 
models. 

In this respect, the model (1.1) is one step closer to the nonlinear NS equation, where the inertial-range behavior 
of structure functions is believed to be related with the Galilean-invariant (and hence built of the velocity gradients) 
operators, which mix heavily in renormalization; see [ |T6| , pT| and references therein. 

Another important question recently addressed is the effects of large-scale anisotropy on the inertial-range statistics; 
see, e.g., Ref. and references therein. In particular, it was shown that in the presence of anisotropic forcing, the 
exponents describing the inertial-range scaling of the passively advected scalar [|4 jsi) and vector [ p3|j29|j30tl fields and 
the turbulent velocity field itself jS^-pSt are organized in a hierarchical order according to their degree of anisotropy, 
with the leading contribution coming from the isotropic sector. The consistency of this picture with the presence of 
nonlocal terms in the equations for the correlation functions, caused by the pressur e co nt ribu tions, has been addressed 
recently and answered positively in Ref. ||2^ on the example of the model (IT)-(O), for the pair correlation 
function in three dimensions, and in Ref. on the example of the general yl-model, for the correlation functions of 
arbitrary order in d dimensions. 

The plan of the paper and the main results are the following. In Sec. |l|, we give the field theoretic formulation 
of the model, its diagrammatic techniques, and derive exact equations for the response function and pair correlation 
function: the so-called Dyson- Wyld equations. 

In Sec. Ill, we study the inertial-range behavi or o f the pair correlation function in the presence of the large- 
scale anisotropic forcing. This issue for the model (L1)~(L3) was already discussed in Ref. |2^, where the numerical 
solutions for the scaling exponents were presented in three dimensions for the isotropic sector and low-order anisotropic 
sectors. 

In this paper, starting from the Dyson-Wyld equations, we give the general recipe of deriving nonperturbative exact 
equations and obtain explicitly transcendental equations for the scaling exponents, related to different irreducible 
representations, in d dimensions. This allows us to give general description of the behavior of the full set of solutions 
in isotropic and anisotropic sectors, and to derive analytical results in all sectors to order 0{e) in d dimensions. Some 
details are given in Appendix A. These results are illustrated by a few nonperturbative solutio ns obt aine d num erically 
in two and three dimensions for the isotropic and low-order anisotropic sectors; later, in Sees. VI D and VIE they are 
confirmed using the RG and OPE techniques and extended to higher-order structure functions. 

In Sec. [V, we perform the ultraviolet (UV) renormalization of the model and derive the corresponding RG equations 
with exactly known RG functions (/? functions and anomalous dimensions of the basis fields and parameters). For 
d^ > 3, these equations possess an IR stable fixed point, which establishes the existence of IR scaling with exactly 
known scaling dimensions of the basis fields and parameters of the model. 

In Sec. we discuss the operator product expansion and its relationship to the issue of anomalous scaling. We 
show that nontrivial exponents describing the inertial-range behavior of the 2n th order even structure function are 
related to critical dimensions of scalar composite operators built of_2n derivatives of the advected field. Explicit 
calculation of the dimensions of relevant operators is given in Sec. 



VI. The case n 



1 can be treated exactly using 
certain functional Schwinger equation (in the case at hand, it has the meaning of the energy balance equation). Like 
in the scalar case the function 5*2 oc r^~^ appears nonanomalous with the simple dimensional exponent; Sec. VIA. 
The case n = 2 is studied in detail; Sec. VI B. The critical dimensions of the relevant family of operators are calculated 
to order 0(e) in d dimensions; they include a negative dimension, and the function ^4 shows anomalous scaling. The 
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families of the anomalous exponents related to the higher-order functions S2n are calculated in Sec. VI C , in the limit 
of large d, for n as high as 6; some technical details are given in Appendi x B. 



Generalization to the case of anisotropic sectors is discussed in Sees. VI D and VIE. There, the RG and OPE 



techniques confirm the general picture established earlier for the pair correlation function (infinite sets of exponents, 
hierarchy, absence of saturation) and extend it to the case of the higher-order structure functions. 

The results obtained are reviewed and discussed in the Conclusion, where the lessons one can learn regarding the 
stirred NS equation and possible generalization to this nonlinear problem are also briefly outlined. 



II. FIELD THEORETIC FORMULATION AND THE DYSON-WYLD EQUATIONS 



The stochastic problem (1.1)-(1.3) is equivalent to the field theoretic model of the extended set of three fields 
$ = {6'',6',v} with action functional 



= 9'Dg6'/2 + 9' [-Vt + voA] 9 - vL>- V/2. 



(2.1) 



The first three terms in Eg. (2.1) represent the Martin-Siggia-Rose-type action for the stochastic problem (1.1), ( |l.2| ) 
at fixed v (see, e.g., ||l^,|l3l and references therein), whi le th e last term represents the Gaussian averaging over v. 



Here Dg and Dy are the correlation functions (1.2) and (1.3), respectively, 9' is an auxiliary transverse vector field, 
the required integrations over x = (t, x) and summations over the vector indices are implied, for example, 



9'dt9 = J dtdx9'^{t,x)dt9,{t,x). 



The pressure term can be omitted in the functional (2.1) owing to the transversality of the auxiliary field: 

dx9'AV = - I dxVdM = 0. 



Of course, this does not mean that the pressure contribution can simply be neglected: the f ield 9' acts as the transverse 
projector and select s th e transverse part of the expression in the square brackets in Eq. (^.l[) 



The formulation (2T) means that statistical averages of random quantities in stochastic problem (L1)-(L3) can be 
represented as functional averages with the weight exp5($), so that the generating functionals of total [G(A)] and 
connected [VF(yl)] Green functions of the problem are given by the functional integral 



G{A) = exp W(A) = j P$exp[5($) -I- 



(2.2) 



with arbitrary sources A = , A^ , A^ in the linear form 

A^^ [dx [Af{x)9',{x) + At{x)9,{x) + AJ{x)v,{x)]. 



The model (2.1) corresponds to a standard Feynman diagrammatic technique with the triple vertex —9'{'vd)9 
9[0jVkVijk with vertex factor 



(2.3) 



where p is the momentum flowing into the vertex via the field 9' . The bare propagators in the frequency-momentum 
(w-p) representation have the forms 



{9,{oj,p)9'^{-u;,-p))^ = {9[{co,p)9,{-Lo,-p))l 



{-iu! + i^op'^)' 



(2.4a) 



{9,{cj,p)9,{-cj,-p))^ = 



(C.2 + ^2p4) 



{9',{u;,p)9'^{-u;,-p))^ = 0, 



(2.4b) 
(2.4c) 
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wher e C ajp) is the Fourier transform of the function Cy(r) from ( |l.2| ); the bare propagator {viVj)o is given by 



Eq. (yD 



The action functional (2.1) is invariant with respect to the dilatation 6 \d, 9' 0' /X, C X^C, where C is the 
correlation function (1^). It then follows that any total or connected Green funct ion w ith n fields 9 and p fields 9' 
is proportional to C'^"^^-'/^. Since the function C appears in the bare propagator (2.4b) only in the numerators, the 
difference n — p is an even non-negative integer for any nonvanishing function; the Green functions with ri — p < 
vanish identically. On the contrary, the 1-irreducible function {9{xi) ■ ■ ■ 9{xn) 9'{yi) ■ ■ ■ 9'{yp))i-ir contains the factor 
(7(p-")/2 g^j^(j therefore vanishes for n — p > 0; this fact will be relevant in the analysis of the renormalizability of the 



model (see Sec. EV). 

The pair correlation functions ($$) of the multicomponent field <& satisfy standard Dyson equation, which in the 
component notation reduces to the system of two nontrivial equations for the exact correlation function I?ij(ti;,p) = 
{9i{uj,p)9j{—Lj,—p)) and the exact response function Qij{Lj,p) = {9i{tJ,p)9j{~uj, —p)) . We shall see below that 
the latter function does not include the correlation function ( |l.2[ ), therefore it is isotropic and can be written as 
Qij(u!,p) = Pij{p)Q{LL>,p) with certain isotropic scalar function 5(w,p). Thus the component equations, usually 
referred to as the Dyson- Wy Id equations, in our model take on the form (cf. Refs. |^ for the scalar and for the 
magnetic models) 



g-\io,p)p,,ip) = [- 



I'aP 



']P.,(p)-Ef/(^,p), 



(2.5a) 



v,,{u;,p) = \g{uj,p)\' a,{p) + ^!',''{io,p) 



(2.5b) 



where ^ and * are self-energy operators represented by the corr espon ding 1-irreducible diagrams; the other 
functions E** vanish identically. It is also convenient to contract Eq. (2.5a) with the projector Pij{p) in order to 
obtain the scalar equation: 



g ^(w,p) = -iw + 



^''iu;,p), 



where we have written 



S'"'(c.,p)^Ef/(c.,p)P,,(p)/(d-l). 



(2.6a) 



(2.6b) 



The feature characteristic of the rapid-change models like is that all the skeleton multiloop diagrams entering 
into the self-energy operators ^ and ^ contain effectively closed circuits of retarded propagators {69') o and 
therefore vanish; it is crucial here that the velocity propagator ( |l.3[) co ntains the i5 function in time and the bare 
propagator ( p. 4c ) vanishes. Therefore the self-energy operators in ( |2.5D are given by the one- loop approximations 
exactly and have the forms 



(2.7a) 



O 



(2.7b) 



The thick solid lines in the diagrams denote the exact propagators {99') and {99); the ends with a slash correspond 
to the field 9', and the ends without a sl ash correspond to 9; the wavy lines denote the velocity propa gato r (|l.3|); 
the vertices correspond to the factor (2^). The first equation does not include the correlation f unct ion (1.2), which 
justifies the isotropic form of the function Qij. The analytic expressions for the diagrams in Eq. (2.7) have the forms 



(d-1) 



^y^y Pky,,3,,(p)p,3,,(p-k)g(c.',p-k) 



(2.8a) 



Ef7(c^,p) = J Vuj' j V\^Vu,^Ap)V,,,,{u',p-'k) 



J»422( P)- 



(2.8b) 



Here we have denoted Vlu' = (iijj'/(27r), used the explicit form ( pT^) of the velocity covariance and the relation 
Pij^i^Ck.) Vi^ji^ip — k) = Piyi^Ck) Vi4ji2(p) for the vertex factor in Eq. (|2.3|). We also recall that the integrations over 
k should be cut off from below at k = ni. 
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The integrations with respect to w' on the right-hand sides of Eqs. ( |2.8| ) give the equal-time response function 
G{k) = J Voj' Q {lo' , k) and the equal-time pair correlation function 2?ij(k) = J 'Doj''Dij{uj' ^'k)\ note that bot h th e 



self-energy operators appear independent of lj. The only contribution to Q comes from the bare prop agator (2.4a) 
which in the t representation is discontinuous at coincident times. Since the correlation function (1.2), which enters 
into the one- loop diagram for ^ , is symmetric in t and i', the response function must be defined &i t — t' by half 
the sum of the limi ts, wh ich is equivalent to the convention G{k) = j Vu' {—ioj' + v^k"^)^^ = 1/2. This allows one to 
write the equation (2.6a) in the form 



where the p-dependent effective "eddy diffusivity" is given by 



(2.9) 



1 - 



(pk)2 



p^k^ - (pk)2 



(d-l)lp-kP 



(2.10) 



It follows from Eq. ( 2.10| ) that the eddy diffusivity can be written as the sum of two parts: ^effip) ~ ^loc~^^non (p), 
where the local part is p-independent and coincides with the expression for the effective diffusivity known in the scalar 
and magnetic cases, while the nonlocal contribution has a finite limit at m = but retains a nontrivial dependence 
on the momentum: 



Do 
2 



Vk 



vq + DoCdm 



■ jd-l) 
2de ' 



(2.11a) 



vnon {p) 



Do 
2 



Vk 



(pk)2 



(pk)^ - p^k^ 
(d-l)p2|p_k|- 



(2.11b) 



Here and below Ca = Sd/{2T:Y and = 2'k'^I'^ /T{d/2) is the surface area of the unit sphere in d-dimensional space 
and r(z) is the Euler Gamma function. The parameter m in z/^^g has arisen from the lower limit in the integral over 
k. For m = 0, equation (2.11b) gives 



Vnon{p) = -DoP V/(47r) 



d/2 



where 



J = 



(c^ + i)r(|)r 



r 1 



sr 1 



r 2 



(2.12) 



(2.13) 



while for p = Q one obtains 



-CdDoni ' 



2d{d + 2)e' 



1^0 + CdDo m~ 



2d{d + 2)e' 



(2.14) 



Equations (2.9)-( ^T4 ) give the explicit exact expression for the response funct ion in our model; it will be used in 
Sec. IV for the exact calculation of the RG functions. The integration of Eq. (2.5b) over the frequency lo gives a 
closed equation for the equal-time correlation function; it is important here that the lo dependence of the right-hand 
side is contained only in the prefactor |tj(w,p)p: 



2v,ff{p)p'V.,{Y>) = Cy(p) + S^/ (p). 
Using Eqs. ( 2.8b| ), ( |2.10D and (2.11), the equation for Vij can be rewritten in the form 



(2.15) 



2vop' V,j (p) ^ C^j (p) ^p'Do j 



Vk 



(pk)2 

p2^2 



Pu^ {p)Vu., (p - k)Pj,j (p) - V^j (p) \ - 2vnon {p)p^ T^.^ (p 



(2.16) 
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The subtracted term in the curly brackets is th e contribution of vi^^ in Eq. (2.15) written in its integral form 
(^.lla|) . For < e < 2, the IR cutoff in Eq. (|2.16| ) can be removed. Indeed, owing to the subtraction, the integral 
over k is finite for m — > 0: the possib le IR divergence at k = is suppressed by the expression in the curly brackets. 
In what follows we set m = in Eq. ( 2.16|). 

It is instructive to compare expression (2.14) for the effective diffusivity in our model with its analogs for the scalar 
and magnetic case s. For those, the nonlocal part of the effective diffusivity vanishes identically, while its local 
part coincides with Eq. (2.11a). Therefore the ratio viod^eff — (d — ^){d + 2)/((i^ — 3) [we have put p ~ and 
neglected the bare diffusivity vq] can be considered as a measure of the nonlocality contribution into the turbulent 
diffusion. It tends to unity as d — > oo, increases monotonically as d decreases and diverges at = 3. This means 
that the nonlocality contribution is negligible at large d (see also Ref. |2^ for the general vector model), becomes 
comparable with the local contribution as d is reduced, and dominates the diffusion in low dimensions (in particular, 
viocl^eff = 5/3 and 4 for d — i and 2, respectively). 



We also notice that, according to Eq. (2.14), the effective diffusivity ^effip) becomes negative for small p ^ m and 
< 3. Therefore the response function in the time-momentum representation, 

where 8(t) is the step function, grows with time for small p, thus signalling that the steady-state solution cannot 
be stable for < 3: any small pert urba tion would lead to the exponential in-time growth of the mean field (9). 
Indeed, one can easily see that Eq. ( 2.15| ) for the correlation function has no solutions at small p and d'^ < 3: its 
left-hand side is negative, while the righ t-hand side is strictly positive for all d. We shall see below in Sec. Ill A that 
the inertial-range solutions of Eq. (2.15) also become singular and disappear in the limit c?^ — > 3. 

Although the instability occurs for unphysical value of d, it deserves some attention as a result of the competition 
between the local and nonlocal contributions: from Eqs. ( ^.11 ) it follows that viq^, is strictly positive for all d, while 
Vnon is strictly negative. In a few respects, such an instability differs from that established in Ref. p7| for the 
magne tic (local) case, in three dimensions, where Vnon = 0, the effective diffusivity coincides with its local part 
( ^.lla ) and is always positive. Thus the response function, and hence the mean (0), show no hint of misbehaviour at 
the threshold of the instability; the latter reveals itself on the level of the pair correlation function and can be related 
to the complexification of the inertial-range exponents; see Refs. [p7-30 . 



III. INERTIAL-RANGE BEHAVIOR AND SCALING EXPONENTS FOR THE PAIR CORRELATION 
FUNCTION IN THE PRESENCE OF LARGE-SCALE ANISOTROPY 



It i s we ll known M-k| that nontrivial inertial-range exponents are determined by zero modes, i.e., the solutions of 



Eq. ( 2.16 ) neglecting both the forcing [C(r) = 0] and the dissipation (i/q — 0). Whatever be the forcing, equations 
for zero modes are linear and SO{d) covariant, and their solution can be sought in the form of decomposition in 
irreducible representations of the rotation group. Equation ( 2.16| ) then falls into independent equations for the 



coefficient functions. In three dimensions, such decompositions were used in Refs. ||2£ 

Below we use more elementary derivation, which allows one to obtain explicitly transcendental equations for the 
scaling exponents, related to different irreducible representations, in d dimensions. We restrict ourselves with the case 
of uniaxial anisotropy, specified by an unit vector n, which is sufficient to find all independent exponents, and use 
explicit expressions for the basis functions in the momentum representation. Then the transversality condition can 
be easily imposed from the very beginning, and there is no need to check it a posteriori. 

In contrast with the real-space Legendre decomposition, employed in Refs. [g^,^, our representation is consistent 
with the rotational symmetry: it can be considered as the projection of the complete SO{d) decomposition onto 
the subspace of the functions with uniaxial symmetry. Therefore, no additional assumptions, like the hierarchy of 
exponents, are needed to disentangle the equations for different anisotropic sectors. 

We start with the isotropic case and then discuss the general situation. 

A. Solution in the isotropic shell 

In the isotropic case, the inertial-range solution is sought in the form 

A,(p) = ^P.,(p)p-'^-'^. (3.1) 



The zero-mode analog of Eq. (2. IE) takes on the form 
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Dn 



1 - 



(pk)2 



|p - kl'^+T p' 



4+1 



fc2 [1- (pk)Vp2p] ■ 

(d- l)|p-k|''+-r+2 



— d — 'J 



(3.2) 



In the following, we shall need the standard reference integrals 



/ 



[i-(pk)Vp2p]" 



(3.3a) 



where 



a + P + d 



d-1 



a + d 



r 



(3.3b) 



The integral ( [3.3a ) is finite in the region of parameters specified by the inequalities a < (convergence at k — s- 0), 



(3 <2n (convergence at k— p ^ 0, improved by the factor [1 — (pk) /p fc 



2L2ln 



\2n\ 



and a + 13 > —d (convergence 



at k — > oo). However, expressio n (|3.3h| ) is meaningful in a wider range of parameters and, in the spirit of analytical 
and dimensional regularizations ||37l , |38| , it can be considered as the analytical continuation of the integral (3.3a) from 
the region in which it converges. The precise meaning of such a continuation is that Eq. (3.3) gives the correct value of 
the integral with proper subtractions which ensure its convergence. For example, if the factor l/|p — kj'*"'"'' is replaced 
with the difference l/|p — k|''+'^ — (that is, th e ze roth term of the Taylor expansion in k is subtracted), the 

integral becomes convergent for a < 2 and expression ( |3.3| ) gives the correct answer for this "improved" integral. 

One can easily see that Eq. ( [3.2[ ) invo lves suc h a subtraction, which improves its convergence at small k. The refore, 
one can use the formal expression (3^) in ( ^^ ) and simultaneously omit the subtracted term. Then Eq. (3_^) takes 
on the form 



/i(£,7)-/2(e-2,7 + 2)/(d-l) 



-2J, 



(3.4) 



with J from Eq. ( 2.13 ). We sh all s ee below that the leading admissible solution of this equation is 7 = 2 — e, so that 
the integrals entering into Eq. (3.2) are convergent for all < e < 2 and the above procedure is internally consistent. 
In Eq. (3.4), we omit the overall nonvanishing factor 



r 1 



l)r('-|ir 



-f + d + e 



e(d+l)r 1 



r 1 



d + 7 



r 2- 



7 + e 



2 / V 2 y V 2 

and obtain the desired equation for the zero-mode exponents in the isotropic case: 



r 1 



(d-l)(d + 7)(2-7-e)/(d+l)+e 



ri 



r 1 



d + 7 



r 2- 



e + 7 



r 2 



r 1 



d + £ + 7 



(3.5) 



The transcendental equation ( [3.5| ) has infinitely many solutions. This means that the inertial-range behavior of the 
correlation function is given by an infinite sum of power like contributions of the form (3.1); the leading term is given 
by the minimal exponent 7. 

Some solutions can be ruled out as not admissible admissible solutions are non- negative for e = (see, e.g., 

[^,Q). The remaining solutions, all having the forms 7 = —d— 2k + 0{e) with non-negative integer fc, are also 
meaningful and describe the behavior of the correlation function at large scales r ^ L [Q ; we shall not discuss them 
in the following. 



It then follows from (3.5) that the leading admissible inertial-range solution is 7 = 2 — e (no corrections of order 
and higher). In the coordinate representation, this corresponds to S2 oc r^^^, that is, the second-order structure 
function is nonanomalous like in the scalar model |^ . 

For small e, all the sublcading exponents can be written in the form 



Ik 



2k 



id-l){d + 2) 
id^-3) 



id+l){d + 2) 



2{d? - 3)2 



1)/C 



{d+1) 



{k^l){d + 2k) 



(3.6) 
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with fc = 2,3,4 and so on. The functions ipiz) = d\'aT(z)/dz can be ehminated from the coefficient 



/C = ^(/c + d/2) - V'(2 + d/2) + V'(fc) - V'(l) 
using the relation -0(2 + 1) = ip{z) + l/z. In order to pre vent the appearance of ambiguities in the zeroth order in s, 



it is convenient to rewrite the right-hand side of Eq. (3.5) such that the Gamma functions have no poles at jk — '^k. 
This is easily done using the relation r(l + 2;)r(l — z) — 7rz/sin(7rz) and is also useful for the numerical solution. 

Nonperturbative solutions of Eq. (3.5) can only be obtained numeric ally . They are illustrated by Fi g. H for d = 2 
(left) and d = 3 (right); the latter is in agreement with Fig. 2 from Ref. g^. One can see that solutions (|3.^ ) exist for 
all < e < 2, decrease monotonically as e grows and turn to 7fc = 2fc — 2 at e = 2. The exponent 7 = corresponds 
to the solution 6ijS{p) which exists for all d (see below). 

It is easily seen that the integrals entering into Eq. (^^) are divergent on these solutions at k — p ^ 0. One 
can say that the knowledge of the exponents is insufficient to discuss the convergence: it is necessary to know the 
behavior of the entire solution in th e re gion of small momenta, |k — p| <C m, where it no longer reduces to a sum 



of power terms with the exponents (3^). However, it is intuitively clear that the form of the solution at such small 
momenta is irrelevant for the calculation of the inertial-range exponents. Indeed, we made no assumptions about the 
form of the solution in that range, sought them in purely powerlike form and managed to derive closed equations for 
the exponents using the prescriptions of the analytical regularization. A simple justification of this procedure follows. 

In coordinate representation, the solution is sought in the form P(r) cx r''C{mr) (for simplicity, here and below we 
omit its vector indices). The convergence problems arise if 7 > 0, which is implied in what follows. It is natural to 
assume that the scaling function C(mr) is such that I?(r) vanishes at r = along with all its derivatives up to the 
71 th order, where n is the maximal integer satisfying the inequality 7 > n. This gives the set of integral relations 

rfqfti • ■ -ftfc 25(q) = 0, A: = 0, 1, . . . ,71 (3.7) 



for the correlation function I?(q) on q '^C{q/m) in momentum representation. The integral entering into Eq. ( |3.2| ) 
can symbolically be written as 

' dqF{q,p)V{q), (3.8) 

where we have introduced the new variable q = p — k. The form of the kernel F(q, p) is clear from the comparison 
with Eq. (t3.2[); the divergence can aris e fro m the region q ^ 0, where the solution behaves as I?(q) c>c q^'^^''. Owing 



to relations ( |3.7j ), the value of integral (|3.8| ) does not change if one subtracts from F(q, p) the first terms of its Taylor 
expansion up to the ?ith order: 

(ncp, - no,P) - - ... - 1, -'^f^] ^(0). (3-8) 

L OQi nl dqi^ ■ ■ ■ dgj„ J 

Now one can set tti = in the function I?(q), that is, replace the exact solution with its inertial-range asymptotic 
expression I?(q) oc q^'^ '': the possible divergence at q = is suppressed by the expression in the curly brackets, which 
behaves as ^"■+2 for q ^ 0. Therefore, one can use the formal rules of analytical regularization and simultaneously 
omit the subtracted terms p8| : this gives the correct answer for the convergent integral with proper subtractions in 



Eq. (3.9) 



We thus conclude that the exponents ( p.q ) may appear in the full solution as correction terms; in the corresponding 
integrals exact solution can be replaced with its powerlike asymptote and the resulting integrals are properly given 
by the rules of analytical regularization. 

Our conclusions are in agreement with those drawn in Ref. |2^ for the equation in coordinate space, although 
the analysis in momentum space appears rather different. In particular, the momentum-space analysis reveals the 
close resemblance between the scalar and vector models: for the former, the correlation function in momentum 
representation also satisfies an integral equation and the above discussion is needed to fix the convergence problem. It 
is also worth noting that the procedure employed in Ref. for the calculation of divergent integrals involves analytical 
continuation from the region of convergence, a nd i s therefore close to the concept of analytical regularization. 



Furthermore, the "realizability" of solutions (3.6) is guaranteed by the fact that in the RG approach they are identi- 
fied with the critical dimensions of composite operators entering into the corresponding operator product expansions; 
see Sec. 0. 



Besides powerlike solutions, Eq. (3.2) also possesses the solution of the form 5ij5{\)). To demonstrate this, we use 



the well-known representation of the d-dimensional delta function 
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S{p) = lim y VxiAx)-" exp[i(px)] = p'"^ \iin[a{p/Ay] 



(3.10) 



and substitute it into Eq. rt3.3). For small a, the integral on its left-hand side is finite. Therefore it vanishes as 



CT — > in Eq. (3.10), and the equation is satisfied. Another such solution, [Sij — dniUj] 6{p), belongs to the first 



anisotropic sector. The both solutions are automatically transverse owing to the relation PiS{p) = 0. In coordinate 
representation, they correspond to constant terms, so that the exponent 7 = can formally be assig ned to them. 
They are indeed present in the pair correlation function, but disappear from the structure function (|_^) owing to the 
relation 5*2 = for r = 0. 

We have already established in the end of Sec. || that the model (|l.l[) -(pT^) becomes unstable at d = \/3 and no 
steady-state solu tion for the correlation function exists in the range of small momenta, p <^ m. It is easily seen that 
expressions (3.6) for the inertial-range exponents diverge in this limit. The same singularity occurs in the first-order 
expressions for another exponents in our model; see, e.g., Eqs. (3.21)-( 3.25| ), ( |6.25| ) and ( |6.33 ) in subsequent sections. 
More over , the RG analysis shows that the actual expansion parameter is 5, oc s/{(P — 3) rather than e itself; see 
Eq. (|4.9|) in Sec. IV. Therefore, the higher-order terms of the e expansions become more and more singular as d 



approaches \/3 from above and any finite-order approximation cannot be trusted. 

Numerical solution shows that for d close to the threshold, the behaviour of exponents 7^ consists of two pronounced 
stages. At the beginning, the exponents decrease rapidly as e increases in agreement with the first-order expression 
in Eq. ( p.6[ ) . Then, for a very small value of e (which tends to zero as d ^ \/3) the instantaneous stabilization takes 
place at an almost constant value 7^ ~ 2fc — 2, and for e = 2 one has 7^ = 2A: — 2 exactly. 

For d — \/3 the solutions do not exist. This fact can be naturally explained, and extended to the other exponents, 
on the basis of the e expansion. The series in e can be rewritten as series in oc £/(rf^ ~3), with coefficients regular at 
d^ = 3. Then they can be inverted into expansions of in 7 [more precisely, in 7 — 2/c = 0(e)], where the coefficients 
are also regular. In other words, the inverted e series can be written as e = [d^ — 3)/(7), where / is a function regular 
as d ^ Vs. For d = \/S this gives e = 0, that is, no solution exists for 7 if e > 0. 



B. Scaling behavior in anisotropic sectors 

Now let us turn to the case of uniaxial anisotropy, specified by an unit vector n in the correlation function ( |1.2| ), 
and introduce the following irreducible Ith rank tensors 

^X7^7'K•••T^,J. (3.11) 

Here and below, XTZV denotes the irreducible part, obtained by subtracting the appropriate expression involving the 
Kronecker delta symbols, such that the resulting tensor is traceless with respect to any pair of indices. In particular, 

Afj;^^ = Ui, Af-p = Uirij - Sij/d, Af- fl = riinjUk - (SijUk + SikUj + Sj^ni)/ (d + 2) and so on. 
It is easy to see that structures ( 3.11 ) are orthogonal on the sphere: 

dnA/'(')AA(") = for I ^ s, (3.12) 

where dn is the area element of t he u nit sphere in the d-dimensional space. 
The contraction of the tensor (|3.1l| ) with a fixed vector p gives 

P^^ ■■■Ph^p' Pliz). z = (np)/p, (3.13) 

where Pi{z) = z' -I- 0(z'~^) are the Gegenbauer polynomials, which reduce to the Legendre polynomials and trigono- 
metrical functions in three and two dimensions, respectively p^ . 

The relation ( 3.12| ) implies that these polynomials arc orthogonal on the sphere even if their arguments are different: 

j dnPi{z)Ps{z') = Q for I ^ s, (3.14) 

where z = (np)/p and z' = {n p')/p ' with any fixed vectors p and p'. 

In terms of the structures ( |3.11 ) and transverse projectors Pij the desired decomposition of the pair correlation 
function can be written as follows: 
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A,(p)=g2?«(p), 



(3.15) 



1=0 



where the summation runs over all even values of I and the coefficient functions have the forms 



Whatever be the coefficient functions A^''^ and B^''\ dependent only on p = |p|, the expression ( [3.16 ) is symmetric in 
the tensor indices i, j and orthogonal to the vector p: piV^''\p) = 



The first tensor structure in Eq. ( 3.16| ) can be expressed in terms of the Gegenbauer polynomials using the relation 
(3.13), while the second structure can be expressed in terms of the polynomials Pi{z) and their derivatives using the 
relation 



l(l-^)^flX^,^,■■■^,P^^P^.■ 



92 



p'Piiz) 



which is obvious from the definition (3.11) and the the relation (3.13). 

Substituting the series ( |3.15D into the zero-mode analog of Eq. ( 2.16| ) then gives the equation for the coefficient 
functions B^''\ whi ch ca n be symbolically written as Ay(p) = 0. Its left-hand side can be decomposed using the 
same tensor structures ( 3.1l| ), so that the equation can be reduced to an infinite family of the scalar equations 



j dnAi,(p) (pn)' = 0, j dnA^ (p)ninj (pn)' ^ = 0, 



(3.17) 



with the integrati on ov er the unit sphere; cf. ( 3.12 ). Owing to the generalized orthogonality relation ( 3.14 ), the /th 
pair of equations (3.17) involves only the coefficient functions A'^^\ i?^'^ with the same index Z, so that the equations 
for different values of / have decoupled. 

In the inertial range, the solution in the I th sector is sought in the form ') = aip-'^-^-'f' , B(') = bip-'^-^-f\ which 
corr esponds to I?*^'-'(r) (x r^' in coordinate representation; for the isotropic shell, / — 0, this reduces to the single term 
(3.1). For I >2, one obtains a pair of linear equations for each pair of coefficients ai, bf. 



0, 



m + oi 



0, 



(3.18) 



where the coefficients C^'^ depend on e, d and the imknown exponent 7;. They all can be expressed in terms of the 
basic integrals (3^). In particular, for / = 2 one obtains 



Cf^ = {d-lfh-{d^ -l)h+dh, 
C''^ =C^2i =-{d-l)h + {d+\)l2 
Cif = id- l)(d2 - 2)7i/2 - {d^ - l)h + dh 



dh 



(3.19) 



where h = Ii + 2J with J from (UJn) and /„ = /„(72 - 2,e + 2). The Ith pair of equations involves integrals 

(I) 

1,1 = Inili — 2, e + 2) with n as high as 1/2 -f 2. The coefficients C^^ for higher values of ^ up to / = 12 are given in 

Appendix A. 

The desired equation for 7^ is obtained as the requirement that the linear homogeneous system (3.18) have nontrivial 
solutions: 



detlC 



(0| 

a/3 I 



0. 



(3.20) 



For any given /, the determinant can easily be written down in terms of standard integrals, provided the coefficients 
are known [see Eq. ( 3.19| ) for / = 2 and Appendi x A for Z < 12]; then using the explicit expressions (3^) for /„ one 
obtains transcendental equations similar to Eq. ( [}.5[ ) but much more cumbersome. We shall not write them down 
explicitly for the sake of brevity and turn to the corresponding solutions. 

It is also worth noting that for Z > 4 the formal divergen ce of t he integrals in Eq. ( |3.2D occurs already for the leading 
exponents and the discussion, similar to that given in Sec. Ill A, is needed to justify the use of analytic regularization. 
For d — 2, the two structures in Eq. ( 3.16 ) coincide and the determinant det |C^'^| vanishes identically. We shall 



return to this case in the end of the Section, and for now on we assume d ^ 2. 
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For any given I > 2 and small e, all possible solutions for the exponents 7; can be written in the form 7; = 
{1 — 2 + 2k) + 0{e) with fc = 0, 1, 2, . . .. The leading exponent corresponds to fc = 0; it is unique for any I > 2 and has 
the form 



7, = (/ - 2) 1 + e 



{d + 2)il- 3){d^ ~4d + 2ld + P - 5/ + 4) 
((i2 -■i){d + 2l- 6)id + 21 - A)id + 21 - 2) 

For all A; > 1, there are exactly two solutions. For fc = 1, they can be written as 

(d + 2) x± 



where the slopes x± satisfy the quadratic equation 



(3.21) 



(3.22) 



v^{d + 21 ~ 4)(d + 21) - [d'^ + (4/ - b)d^ + A{f - 4? + 2)d^ + (-4?^ - 2f + 14/ - A)d - 2l{l - 1){1 - 2)(3Z + l)]x - 

-lil - l)(d + / - 1) [d^ + (3/ - 5)^2 + (2/2 - 11/ + 8)d + 2(-/^ + /^ + 3/ - 2)] = 0. 

In particular, for / = 2 this gives 

x± = {d^ + - 8d - 16 ± ^/{d + 4)(d5 + 2d^ - 7d^ - M"^ + 8d+ 16)} / 2{d + 4). (3.23) 
For fc = 2, the solutions have the form 



7, = (/ + 2) + e 



(d + 2)[-d-^ + (2 - 6l)d^ + (3 + 9/ - 11/2)^2 + (-6 + 6/ + 16/^ - 6P) d + 1{P - 1){1 + 10)] 



(d2 - 3)(d + 2/ - 2)((i + 2/)(d + 2/ + 2) 



(d-2)(d + 2) 2 
7; = (/ + 2) - e — _ — + 0(e ). 



The situation simplifies for all fc > 3: then the solutions in order 0{e) are degenerate and have the form 

7. = (i-2 + 2.)-,. '-'-/-' + ^) +0(e'). 



(3.24) 



(3.25) 



that is, the slope is the same as for the second solution in Eq. ( ^.24 ) and the stand ard s lope ( ]3.6| ) for the isotropic 
sector. However, the degeneracy is removed by the 0{e'^) terms, and the two solutions (3.25) do not coincide identically. 

At the opposite edge, s = 2, all the solutions can also be found analytically for any given / > 2. They take on the 
values / — 2 (single), /, 1 + 2, 1 + 4, and so on (two-fold degeneracy for d ^ 2 and single otherwise). In addition to 
these "standard" values, for all / > 2 there are exactly two d-dependent solutions 7^^; they satisfy certain quadratic 
equations and have the forms 



7.+ = 



1 \d^ + 4/(i2 - d2 _ _|. 4/2^ _ ^ 4^2 

~2\ 



(d-1) 



1 /d3 + 4/d2-9d2-4/d + 4/2^ + 8-8/ + 4/2 
~2\\ 



(d-1) 



(3.26) 



(only one solution of each equation is admissible). Note that 7^^ > 7^^ for all / and c? > 1. For large d, exponents 
( |3.26 ) behave ss = I + 0{\ j d) and — I ~ 2 + 0(\l d), respectively, so that all solutions become "standard" at 
d = 00. 

Nonperturbative solutions for intermediate values of e between and 2 can only be obtained numerically. We have 
performed the calculation in two and three dimensions for / < 12; the results are illustrated by Fig. || for / = 2, 8 and 
12. For d = 3 and / < 10, our solutions are in agreement with the results presented in Fig. 2 of Ref. |^6|, except for 
the case / — 2: the behaviors of the solutions 72 = 2 + 0(e) and 72 = 4 + 0(e) are different. We believe that this 
disagreement is not conceptual and is explained by calculational errors in ||2^ . For the sake of brevity, we do not give 
the solutions for / = 4, 6 and 10, which in three dimensions are in agreement with 1 26 . The exponent 7 = for / — 2 
corresponds to the solution — driiTij] 6{p), which exists for all d (see Sec. Ill A). 

The figures illustrate the following qualitative behavior of the solutions, which holds for all d ^ 2 and / > 2. The 
leading solution (3.21) exists for all e and tur ns to / — 2 for e = 2. In fact, it is hardly distinguishable from a constant, 
7/ w / — 2, for all values of e (72 = 0, see Sec. IIIA). 
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For any I, some critical value kc = kc{l,d) exists such that for all k > kc the behavior of the solutions 7; = 
I ~ 2 + 2k + 0(e) is simple: the both solutions exist for all e, decrease slowly as e grows and turn to 7/ = Z — 4 + 2fc 
for £ = 2. (In fact, the both solutions corresponding to given / and k are very close to each other for all values of e.) 

For fixed I, the critical value kc decreases as d increases, so that all solutions with k = 1,2,... become simple (in 
the above sense) provided d is large enough. For fixed d, the critical value k^ increases with I; in particular, in three 
dimensions kc = 2 for I ~ 2, 4, 6, kc ^ 3 for Z = 8, 10 and kc — 4 for I = 12. 

An interesting interaction between the solutions with a fixed I and different k occurs for 1 < k < kc- Two branches 
starting at e = with different values of k can coalesce and disappear for some value of e between and 2. Another 
possible process is the creation of a pair of solutions for some < e < 2. A solution that starts at e = can 
annihilate with a solution from a pair that was created for some finite value of e. The interplay between these 
creation-annihilation processes can produce a very complicated pattern, as illustrated by Fig. ^for the sectors I = 2, 
8 and 12 (see also Fig. 2 in Ref. |6|). 

It turns out, however, that the creation and annihilation of solutions eventually compensate each other in the sense 
that the number of branches starting at e = with < k < kc is equal to the number of branches arriving at e = 2 
and confined between the leading solution (7; = Z — 2 for e = 2) and the lowest "simple" s olutio n (7; = I — 4 + 2kc 
for e = 2). The balance is possible owing to the existence of two "nonstandard" solutions (3.26) at the edge e — 2. 
They also determine the boundary between the solutions with "simple" and "interesting" behavior: the uppermost 
solution with the interesting behavior, 7; =1 — 2 + 2kc + 0{e), turns to j'j^ at e = 2. [Some reservations are needed if 
a "standard" solution at e = 2 lies between the roots or coincides with one of them. In particular, for c? = 3 and 
1 = 2, the standard solution 72 = 2 exists but is isolated in the sense that no real branches attach it from the region 
e < 2. For d = 3 and I = 12, one obtains 7j^ — 16, and this standard value acquires three-fold degeneracy]. 

The behavior eventually simplifies in the limit d — > 00. All the solutions become simple in the above sense and they 
are described by straight lines: 7; = Z — 2 for the leading solution and 7/ = Z — 2 -t- 2A; — e for all A; > 1. 

The annihilation of coalescing solutions actually means that they become complex: the effect known for the magnetic 
model where it occurs in the isotropic shell. It was argued in Refs. [ p7|p9|] that the complexification leads to 
the instability of the steady state (exponential growth of the pair correlation function). We shall not discuss this 
important issue here and only stress an essential distinction between the two cases. In the magnetic model, the leading 
admissible exponent 7 = 0{e) coalesces with the solution 7 = —d + 0(e) , wh ich is not ad miss ib le a nd describes the 
large-scale behavior at r ^ L [see the remark and references below Eq. (3.5)]. In model (1.1)-(1.3) the coalescence 
occurs only in anisotropic sectors and only for nonleading admissible exponents. If the steady state remains stable, 
the inertial-range behavior in the corresponding sectors will include oscillations on the powerlike background; cf. the 
discussion in Ref. []l3|| . 

In two dimensions, the tensor structures in decomposition ( |3.16| ) become coincident, and the determinant det |C^'^| 

in Eq. ( [3.20| ) vanishes identically. All the coefficients C^'^ for d = 2 become equal up to the sign [see Eq. (p^ for 
1 = 2]. Therefore, the equation for the exponents ji can simply be written as 



'-^11 



0. 



(3.27) 



All solutions of the d-dimensional equation (3.20) have well-defined limits as d — s- 2, and all true two-dimensional 



solutions are indeed recovered in this limit. However, this limit gives more solutions than the correct two-dimensional 
equation: one half of the solutions obtained in the limit d — s- 2 from the d-dimensional case do not satisfy Eq. (3.27) 
and should be discarded. This behavior is illustrated by Fig. ^, where the solid lines on the diagrams with d = 2 and 
I > 2 denote solutions obtained both from the two-dimensional equation ( 3.27| ) and as limits d — > 2 from d-d imens ional 
solutions, and the dashe d line s denote spurious solut ions w hic h are obtained in the limit d ^ 2 from Eq. ( 3.20| ), but 
do not satisfy equation ( 3.27 ). We shall see in Sees. VI C and VID that similar effect is encountered in the RG and 
OPE approach: all critical dimensions have well-defined limits as d — *■ 2, but some of them should be ruled out due 
to linear relations between composite operators which hold in two dimensions. 



IV. RENORMALIZATION, RG FUNCTIONS AND RG EQUATIONS 

The analy sis of the UV divergences is based on the analysis of canonical dimensions |3^,^^. Dynamical models 
of the type ( ^.l[ ), in contrast to static models, have two scales, so that the canonical dimension of some quantity F 
(a field or a parameter in the action functional) is described by two numbers, the momentum dimension dp and the 
frequency dimension d"^. They are determined such that [F] ~ [L] [T] , where L is the length scale and T is 
the time scale. The dimensions are found from the obvious normalization conditions d^ = — d^ = 1, d^ = d^ = 0, 
dfj = dj' = 0, d^ = — d^ = 1, and from the requirement that each term of the action functional be dimensionless (with 
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respect to the momentum and frequency dimensions separately). Then, based on dp and dp, one can introduce the 
total canonical dimension dp = dp + 2dp (in the free theory, dt (x A), which plays in the theory of renormalization 
of dynamical models the same rol e as the conventional (momentum) dimension does in static problems |Q . 

The dimensions for the model (2^) are given in Table |, including the parameters which will be introduced later 
on. From the Table it follows that the model is logarithmic (the coupling constant go is dimensionless) at e = 0, so 
that the UV divergences have the form of the poles in e in the Green functions. 

The total canonical dimension of an arbitrary 1-irreducible Green function F = ($ • • • <l')i-ir is given by the relation 



dr = rfr 



2d'^=d 



2 - N^d. 



(4.1) 



where iV$ — {Ngi , Ng, N^} are the numbers of corresponding fields entering into the function F, and the summation 
over all types of the fields is implied. The total dimension dr is the formal index of the UV divergence. Superficial 
UV divergences, whose elimination requires counterterms, can be present only in those functions F for which dr is a 
non-negative integer. 

Analysis of the divergences should be based on the following auxiliary consid erat ions [i^JTtI] : 

(i) From the explicit form of the vertex and bare propagators in the model (2.1) it follows that Ngi — Ng = 2Nq 
for any 1-irreduciblc Green function, where Nq > is the total number of bare propagators {09)^ entering into the 
function (see Sec. |l|). Therefore, the difference Ngi — Ng is an even non- negative integer for any nonvanishing function. 

(ii) If for some reason a number of external momenta occurs as an overall factor in all the diagrams of a given 
Green function, the real index of divergence dp is smaller than dp by the corresponding number (the Green function 
requires counterterms only if dp is a non-negative integer). In the model ( |2.lD , the derivative d at the vertex 9'(yd)9 
can be moved onto the field 9' by virtue of the transversality of the field v. Therefore, in any 1-irreducible diagram 
it is always possible to move the derivative onto any of the external "tails" 9 or 9' , which decreases the real index of 
divergence: dp = dr — iVg — Ngr. The fields 9, 9' enter into the counterterms only in the form of derivatives 96*, 89' . 

From the dimensions in Table | we find dv = d + 2 ~ N^ + Ng ~ {d + l)Ngi and dr = (d -h 2)(1 -Ng,)- N^. It then 
follows that for any d, superficial divergences can only exist in the 1-irreducible functions {9' 9 . . . 6')i_ir with Ng, = 1 
and arbitrary value of Ng, for which dr = 2, dp = 0. However, all the functions with Ng > Ng, vanish (see above) and 
obviously do not require counterterms. We a re left with the only superficially divergent function {9'9)i^-„, which does 
not depend on the correlation function ( |l.2| ) and therefore is isotropic; see Sec. |l[ The corresponding counterterm 
must contain two symbols d, and owing to the isotropy and transversality conditions reduces to the only structure 
9'M. 

Inclusion of this counterterm is reproduced by the multiplicative renormalization of the parameters g^, in the 
action functional ( ^.l| ) with the only independent renormalization constant Z^,: 



^0 



Z„ — z,, 



(4.2) 



Here is the reference mass in the minimal subtraction (MS) scheme, which is always used in what follows, g and v are 
renormalized an alog s of the bare parameters g^ and vq, and Z = Z{g,e,d) are the renor mal ization constants. Their 
relation in Eq. (|4.2| ) results from the absence of renormalization of the last term in Eq. (2T). No renormalization of 
the fields and the "mass" m is required, i.e., Z^ = 1 for all and mo = m. The renormalized action functional has 
the form 



Sni^) = 9'Dg9'/2 + 9' [-V* + j/Z^A] 9 - vD-\/2, 



(4.3) 



where the amplitude Dq from Eq. (1.3) expressed in renormalized parameters using Eqs. (4.2): Dq = g^VQ = g^J!^v■ 

The explicit form of the constant Z^ is determined by the requirement that the 1-irreducible function {9'9)i^-„ 
expressed in renormalized variables be UV finite (i.e., be finite for e ^ 0). This requirement determines Z^ up to 
an UV finite contribution; the latter is fixed by the choice of the renormalization scheme. In the MS scheme all 
renormalization constants have the form "1 -I- only poles in e." The function {9'9)i--„ in our model is known exactly; 
see Eqs. (p?9|)-(2.14) in Sec. ||. We substitute Eqs. (4.2) into it and choose Zy to cancel the pole in e in the resulting 
expression. This is equivalent to the requirement that Veff{p) be finite; its pole part is independent of p and is 
therefore contained in Eq. ( ^.14| ). This gives: 



Z,, 



(rf^-3) 
2d[d + 2)e' 



(4.4) 



with coefficient Cd from (2.11). The result (4.4) is exact, i.e., it has no corrections of order g , g^, and so on; this is 
a conseque nce o f the fact that the one-loop approximation (2.7a) for the response function is exact. Also note that 
expression (4.4) differs from the exact expression for Z^ in the scalar M and magnetic [p3l cases. 
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The relation Cq) = e, /i) (where eo = {go, t'o, is the complete set of bare parameters, and e = {g, v, m} 

is the set of their renormalized analogs) imphes W{A, eo) = Wr{A, e,fi), where W is the functional (2^) and Wn is 
its renormahzed counterpart obtained by the replacement S Sr. We use I?^ to denote the differential operation 
fidf^ for fixed eo and operate on both sides of this relation with it. This gives the basic RG differential equation: 



(4.5) 



where Vrg is the operation Va expressed in the renormalized variables: 



V 



RG 



V,,+l3{g)dg-^,{g)V,. 



(4.6) 



In Eq. (4.6), we have written 'Dx = xdx for any variable x, and the RG hmctions (the (5 function and the anomalous 
dimension 7) are defined as 



"frig) — In Zp for any Zp, 



(4.7a) 



l3{g) = V^g^g[~e + -iM]. 



(4.7b) 



The relati on betwe en (3 and 7 in Eq. ( 4.7b ) results from the definitions and the last relation in Eq. (12). From the 
relations (4.4) and (4.7) one obtains explicit expressions for the RG functions: 



7.^(5) 



-eVg \nZ^ 
- Vg In Z„ 



gCd 



2d(d + 2)' 



From Eq. ( 4.7b| ) it follows that the RG equations of the model have an IR stable fixed point [/3(g* 
with the coordinate 

_ 2d{d + 2)e 
^* ~ Cd{d^~3)- 



(4.8) 

= 0, P'{g,) > 0] 
(4.9) 



From the relation between the RG functions in Eq. ( 4.7b| ) the value of 7i/(.g) at the fixed point is found exactly: 
it = lu{g*) = e- 

For d^ < 3, the fixed point is negative and therefore not accessible for the RG flow with physical (positive) initial 
data for g. This is in agreement wit h the conclusion of Sec. ^ that no stable steady state exists in the model for 
< 3 (see the discussion below Eq. ( 2.16 )). 

For d^ > 3, the fixed point is positive; this establishes the existence of scaling behavior in the IR region (Ar ^ 1 
and any fixed mr) for all correlation functions of the model. Let F be some multiplicatively renormalized quantity 
(say, a correlation function involving composite operators), i.e., F = ZpFr with certain renormalization constant Zp. 
It satisfies the RG equation of the form \Drg + 1p]Fr^ = with 7i? from ( 4.7a ) and Vrq from (4^). The solution of 
the RG equation then shows that in the IR region F takes on the scaling form 



F ~ A-'^^ r- 



'(,p{mr), 



where 



Ap = d''p + A^d'^+-f*p, A^ = 2-7* 



(4.10) 



(4.11) 



is the critical dimension of the function F, d^ and dp are its frequency and total canonical dimensions, ^p = lp{g*) 
is the value of its anomalous dimension at the fixed point, A^^ = 2 — 7* = 2 — e is the critical dimension of the 
frequency, and £,p(jnr) is the scaling function whose form is not determined by the RG equation itself. Derivation of 
Eq. (4.11) and more detail can be found in Refs. [p^| jl^ j l7|,|3^ . In particular, for the structure functions (1.6) with 
dp = 0, d% ~ —n/2 (see Table |) and jp — (see Sec. VI) one obtains 



5„(r) = D, 



- n-"/2^«(l-e/2) 



(4.12) 



so that the dependence on the UV scale A disappears, while the dependence on the IR scale m is contained in the 
scaling functions ^„(mr). 



15 



V. OPERATOR PRODUCT EXPANSION AND ANOMALOUS SCALING 



Representations ( }4.10D -(4.12) for any scaling functions ^(mr) describe the behavior of the correlation functions 
for Ar ^ 1 and any fixed value of mr. The inertial range corresponds to the additional condition mr <^ 1. The 
form of the functions S^{fnr) is not determined by the RG equations themselves; in analogy with the theory of critical 
phenomena, their behavior for mr — s- is studied using th e op er ator product expansion (OPE); see Refs. p7| , ^ . 
Below we concentrate on the equal-time structure functions (1^), ( 4.12 ). 

According to the OPE, the behavior of the quantities entering into the right-hand side of Eq. ( |l.6| ) for r = x— x' ^ 
and fixed x x' is given by the infinite sum 



(5.1) 



where Cp are coefficients regular in m? and F are all possible renormalized local composite operators allowed by the 
symmetry. More precisely, the operators entering into the OPE are those which appear in the naive Taylor expansion, 
and all the operators that admix to them in renormalization. 

In what follows it is assumed that the expansion is made in irreducible tensors (scalars, vectors and traceless tensors); 
the possible tensor indices of the operators F are contracted with the corresponding indices of the coefficients Cp- 
With no loss of generality, it can also be assumed that the expansion is made in "scaling" operators, i.e., those having 
definite critical dimensions Ap (see Sec. [V]| ). 

The structure functions ( |l.6D are obtained by averaging Eq. (5.1) with the weight expiS^j, the mean values {F) 
appear on the right-hand side. Their asymptotic behavior for m — s- is found from the corresponding RG equations 
and has the form (F) cx m^^ 



From the RG representation ( 4.12 ) and the operator product expansion (5T) we therefore find the following ex- 
pression for the structure function in the inertial range (Ar 3> 1, mr <C 1) 



5„(r) = 



-e/2) 



(mr) (mr) 



A. 



(5.2) 



where the coefficients Ap are regular in (mr)^. 
Some general remarks are now in order. 

O wing to the translational invariance, the operators having the form of total derivatives give no contribution to Eq. 
(5.2): {dF{x)) — d{F{x)) — d x const = (these operators become relevant if the stirring force in Eq. (1.1) violates 
the translational invariance, like in the problem discussed in Ref. [^ ). 

In the model (1.1)-(1.3), the operators with an odd number of fields 9 also have vanishing mean values; their 
contributions vanish along with the odd structure functions themselves (they will be "activated" in the presence of a 
nonzero mixed correlation functi on ( v/) ; we shall not discuss this possibility here) . 



If the tensor Cij(r) in Eq. (O) is taken to be isotropic, the model becomes SO{d) covariant and only the 
contributions of the scalar operators survive in (5.2). Indeed, in the isotropic case the mean value of a tensor operator 
depends only on scalar parameters, its tensor indices can only be those of Kronecker delta symbols. It is impossible, 
however, to construct nonzero irreducible (traceless) tensor solely of the delta symbols. 

In the presence of anisotropy, irreducible tensor operators acquire nonzero mean values and their contributions 
appear on the right-hand side of Eq. (p^). Like in Sec. HI, consider the case of the uniaxial anisotropy, specified by 
an unit vector n in the correlation function (1.2). In this case, the mean value of a ^ th rank traceless operator involves 
the vector n along with the delta symbols and is necessarily proportional to the / th rank symmetric traceless tensor 
from Eq. (3.11). The contraction with the corresponding coefficient Cp gives rise to the Zth order Gegenbauer 
polynomial Pi{z) with z = (rn)/r; see Eq. ( |3.13|) . In general, the expansion in irreducible tensors in (5T) after the 
averaging leads to the SO{d) decomposition employed in Rcfs. p6| , ^ , the I th shell corresponding to the contribution 
of the Zth rank composite operators. 

The feature characteristic of the models describing turbulence is the existence of the so-called "dangerous" composite 
operators with negative critical dimensions; see Refs. [ p"6|JT7| . Their contributions into the OPE give rise to singular 
behavior of the scaling functions for mr 0, that is, the anomalous scaling. The leading term in the I th anisotropic 
sector is given by the Hh rank tensor operator with minimal (not necessarily negative) dimension A[_F]. 

Since A^^ = dp -\- 0(e), see Eq. (4.11), the operators with minimal Ai? are those involving maximu m p os sible 



number of fields 9 and minimum possible number of derivatives (at least for small e). Both the p robl em (1.1)-([L.3|) 
and the quantities ( |l.6|) possess the symmetry 9 9 + const. It then follows that the expansion ( ^.l[ ) involves only 
operators invariant with respect to this shift and therefore built of the gradients of 9. 
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As already mentioned above, the operators entering into the right-hand side of Eq. (5.1) are those which appear in 
the Taylor expansion, and those that admix to them in renormalization. The leading term of the Taylor expansion 
for Sn is the 2nth rank operator which can symbolically be written as (90)"; its decomposition in irreducible tensors 
gives rise to operators of lower rank s. These contributions exist in the OPE (before averaging) even if the stirring 
force in not included into Eq. (1.1); in the language of Refs. |^-|^ it is then tempting to identify them with zero 
modes, i.e., the solutions of the homogeneous (unforced) analogs of the closed exact equations satisfied by the equal- 
time correlations. In the presence of the stirring force, operators of the form [89)'^ with k < n admix to them in 
renormalization and appear in the OPE; their contributions correspond to solutions of the inhomogeneous equations. 
Owing to the linearity of problem ( |l.l| ) , operators with k > n (whose contributions would be more important) do not 
admix in renormalization to the terms of the Taylor expansion for Sn and do not appear in the corresponding OPE. 
All these operators have minimal possible canonical dimension dp ~ (see Table |) and determine the leading terms 
of the mr — )■ behavior in the sectors with j < 2n. Operators involving more derivatives than fields 9 (and thus 
having canonical dimensions dp — 1, 2 and so on) determine correction terms for j < 2n and leading terms for higher 
anisotropic sectors with j > 2n. The renormalization and dimensions of the most important operators are studied in 
the next Section. 



VI. RENORMALIZATION AND CRITICAL DIMENSIONS OF COMPOSITE OPERATORS 

We recall that the term "composite operator" refers to any local (unless stated to be otherwise) monomial or 
polynomial built of primary fields and their derivatives at a single spacetime point x = (t,x); see Refs. p7| , ^ . Since 
the arguments of the fields coincide, correlation functions with such operators contain additional UV divergences, 
which are removed by additional renormalization procedure. For the renormalized correlation functions the RG 
equations are obtained, which describe IR scaling of certain "basis" operators F with definite critical dimensions 
Ap = A.[F]. Due to the renormalization, A[F] does not coincide in general with the naive sum of critical dimensions 
of the fields and derivatives entering into F. As a rule, composite operators "mix" in renormalization, i.e., an UV 
finite renormalized operator F^ has the form F^ = F+ counterterms, where the contribution of the counterterms is 
a linear combination of F itself and, possibly, other unrenormalized operators which "admix" to F. 

Let F = {Fa} be a closed set, all of whose monomials mix only with each other in renormalization. The renormal- 
ization matrix Zp = and the matrix of anomalous dimensions jp = {70/3} for this set are given by 

Fa = Yl = Zp^V^Zp, (6.1) 





and the corresponding matrix of critical dimensions = {Ac,^} is given by Eq. (4.11), in which dp and dp are 
understood as the diagonal matrices of canonical dimensions of the operators in question (with the diagonal elements 
equal to sums of corresponding dimensions of all fields and derivatives constituting F) and 7^^ = lF{g*) is the matrix 



(O) at the fixed point (4.9) 



Critical dimensions of the set F = {Fa} are given by the eigenvalues of the matrix Aj;-. The "basis" operators that 
possess definite critical dimensions have the form 

Fl'^^ ^Y.UapF^ , (6.2) 

where the matrix Up = {Uap} is such that A'^ = UpApUp"^ is diagonal. 

In general, counterterms to a given operator F are determined by all possible 1-irreducible Green functions with 
one operator F and arbitrary number of primary fields, F = {F{x)^{xi) . . . <f>(a;2))i-ir. The total canonical dimension 
(formal index of UV divergence) for such a function is given by 

dr = dp — Af$c?$, (6.3) 

with the summation over all types of fields entering into the function. For superficially divergent diagrams, dr is a 



non-negative integer; cf. Sec. IV. 

Let us begin with the simplest operators of the form 0"(a;), with free tensor indices or involving any contraction. 
From Table | in Sec. ^ and Eq. (|!|) we obtain dp = -n, dr = -n + Ne~ N^- id+ l)Ng,. 

From the analysis of the diagrams it follows that the total number Ng of the fields 9 entering into the 1-irreducible 
function F = {O"^ {x)9(xi) ■ ■ ■ 6'(x„g))i_ir cannot exceed the number of the fields 9 in the operator 6*" itself, i.e., Ng < n 
[cf. item (i) in Sec. L^. Therefore, the divergence can only exist in the functions with A'v — Ngi = 0, and arbitrary 
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value of n = Ng, for which the formal index vanishes, dr = 0. However, at least one of Ng external "tails" of the 
field 9 is attached to a vertex 9'{vd)9 (it is impossible to construct nontrivial, superficially divergent diagram of the 
desired type with all the external tails attached to the vertex F) , at least one derivative d appears as an extra factor 
in the diagram, and, consequently, the real index of divergence dp is necessarily negative. 

This means that the operators 9" require no counterterms at all, i.e., they are in fact UV finite , 9^ = Z [9"]^ with 
Z = 1. It then follows that the critical dimension of 9"'{x) is simply given by the expression ( 4.11 ) with no correction 
from 7^ and is therefore redu ced to the sum of the critical dimensions of the factors: A[9"] — nA[9] — n{—l + e/2). 
Since the structure functions (1.6) are linear combinations of pair correlation function involving the operators 0", this 
relation shows that they indeed satisfy the homogeneous RG equation (l.f;), discussed in Sec. 

In the OPE for the pair correlation function, analogous to Eq. (5.2), the operators 9"^ and 9i9j with the dimensions 
2(— 1 + e/2) give rise to constant terms. They correspond to the solutions with 6{p), discussed in the end of Sec. 
Such terms, caused by various operators of the form 9", are also present in higher-order correlation functions. They 
disappear from the structure functions (1.6), whose inertial-range behavior is determined by operators built only of 
gradients (see Sec. 



A. Scalar operators of the form (89)^ and the scaling of S2 



The leading terms of the inertial-range behavior of the second-order structure function 5*2 are determined by the 
critical dimensions of the composite operators built of two gradients: 



di9jdi9j, 



di9jdji 



(6.4) 

and its dimension 



For the transverse field 9, the second operator reduces to a tot al d erivative, F2 
A2 = 2 + 2/S.0 — e does not appear on the right-hand side of Eq. (5.2). 

The dimension of the first operator is found exactly: Ai = 0. This can be demonstrated using the Schwinger 
equation of the form 



59[{x) 



a,(a;)exp[5fl($)+^$] 1=0 



(6.5) 



(in the general sense of the term, Schwinger equations are any relations stating that any functional integral of a total 
variational derivative is equal to zero , se e, e.g., Refs. |^,^) . H ere Sb. is the renormalized analog of the action (| 



and the notation introduced in Eq. (2.2) is used. Equation (6.5) can be rewritten in the form 



{{9'Dg9 - Vt[0V2] - d,[9,V] + i^Z,A[0V2] - vZ,F^))a = -Ae,5WR{A)/5Ag. 



(6.6) 



Here Dg is the correlation function (|l.2|), 9^ = 9i9i, ((. . .))a denotes the averaging with the weight exp[5fl(<I>) + A<f>], 
Wr is determined by Eq. (2.2) with the replacement S —^ Sr. , and the argument x common to all the quantities is 
omitted. The contribution with the pressure V from Eq. (1.5) arises due to t he fa ct that the differentiation in Eq. 
( |6.5|) is performed with respect to a transverse field; see the remark below Eq. (2.1). 



The quantity {{F))a is the generating functional of the correlation functions with one operator F and any number 
of fields therefore the UV finiteness of the operator F is equivalent to the finiteness of the functional {{F))a- The 
quantity on the right-hand side of Eq. ( |6.6| ) is finite (a derivative of the renormalized functional with respect to a 
finite argument), and so is the operator on the left-hand side. Our operator Fi does not admix in renormalization to 
the operator 9' Dg9 (Fi contains too many fields 9), and to the operators Vt[6'^/2], di[9iP] and A[6'^/2] (they have the 
form of total derivatives, and Fi does not reduce to this form). On the other hand, the operators 9'Dg9 and di[9iP] 
do not admix to Fi (they are nonlocal, and Fi is local), while the derivatives [6*^/2] and A[6'^/2] do not admix to 
Fi owing to the fact that each field 9 enters into the counterterms of the operators Fn only in the form of derivative 
d9 (see above). Therefore, all three types of operators entering into the left-hand side of Eq. (6.6) are independent, 
and they must be UV finite separately. 

Since the operator vZ^Fi is UV finite, it coincides with its finite part, i.e., vZ,jFi — vF-^, which along with the 
relation Fi = ZiF^ gives Zi = Z~^ and therefore 71 = — 7,y. For th e cr itical exponent Ai = e + 7]^ we then obtain 
Ai = exactly (we recall that 7* — e; see the discussion below Eq. (O) in Sec. IV). 

It then follows from (p^) that the leading term of the inertial-range behavior of the second-order structure function 
has the form S2 oc Dq^ r'^"'^, in agreement with the solution 7 = 2 — £ obtained in Sec. Ill from the exact equation. 



Therefore, this function is not anomalous, like its analog for the scalar model ||^-|^, and the anomalous scaling reveals 
itself only on the level of the fourth-order structure function. 
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B. Scalar operators of the form {89)^ and the anomalous scaling of S4 



Let us turn to the scalar composite operators built of four gradients d9, which cannot be reduced to the form of 
total derivatives. This family includes six independent monomials, all of which can be obtained from the fourth rank 
operator ^"^/'^ = diOm djOn dkdp diOg by various contractions of the tensor indices: 



iijk 
jkll ' 



iijk 

jiki ■ 



(6.7) 



At first glance, it seems that one can add another independent monomial, F7 = ^'[jjl^ but in fact it reduces to Fi up 
to total derivatives: 



3Fi - 6F7 = d. 



(6.8) 



where the notation is analogous to that in Eq. (6/7); see Ref. ||4l| . 

Now let us turn to the calculation of the renormalization constants for the family (^^) in the one- loop approximation. 
Let Ta{x ; 6) be the generating functional of the 1-irreducible Green functions with one composite operator from 
Eq. (6.7) and any number of fields 9. Here x = {t, x) is the argument of the operator and 9 is the functional argument, 
the "classical counterpart" of the random field 9. We are interested in the fourth term of the expansion of Ta{x; 9) in 

9, which we denote T'it\x] 9). It has the form 



V'i\x-9)^l 



4! 



dxi 



dXi 9i^{xi) ■ ■ ■9i^{x4) {Fa{x)9i^{xi) ■ ■ ■ 9i^{x4))l-ir- 



In the one-loop approximation this function is represented diagrammatically as follows: 



Tl^)ix;9)^F^ + ^ A . 



(6.9) 



Here the thin solid lines denote the bare propagator {99') q from Eq. ( 2.4a ), the ends with a slash correspond to the 
field 9', and the ends with out a slash correspond to 9] the wavy line denotes the velocity propagator ( |1.3| ); the vertices 
correspond to the factor (2^). The first term is the "tree" approximation, and the black circle with two attached 
lines in the diagram denotes the variational derivative 

V};tli^; XI, X2) = d^F^{x)/69,^ixi)d9,,ix2). (6.10) 
The diagram is written analytically in the form 

J dxi ■ ■ ■ J dxiVlf^{x; Xl,X2){9^,{xl)9l^{x3))o{9,^{x2)9'^^{x4))o{v,.^{x3)v,^{x4))o^,.J,^{x3)^,^9,^{xi), (6.11) 



with the bare propagators from E qs. ( 1.3) and (2.4); the derivatives appear from the ordinary vertices ( p.3| ). It is 
convenient to represent the vertex (3.10) in the form 



^iilii^'^ a^i,a;2) = " '^^"^ di^di^i^{x - xi) diji^j^ix - X2), 

^^isiA^^i^ie 



d^Fa^ia) 



(6.12) 



where 



d,j (x - x') = 6{t - t')P,j (x - x') = S{t - t') j PkPy (k) exppk • (x - x')] 

is the delta function on the transverse subspace. The first (combinatorial) factor in Eq. (6.12) is understood as follows: 
the gradients di9j{x) in the operator F^ are replaced with a constant tensor , the differentiation is performed with 
respect to its elements, and after the differentiation they are replaced back with the gradients, aij — > di9j{x). 

Using the identity dkSij{x — x') — —d'f.5ij{x — x') and the integration by parts, the derivatives can be moved from 
the vertex onto the propagators, and the integrations with respect to xi and X2 are then easily performed: 



^^ijii ^^i%i2 



dx3 / dx4{di^9i^{x)9'^ (x3))o{di^9i^{x)9'i^{x4))o{vir^{x3)vi^{x4))odi^di3{x3)dig9i^{x4). (6.13) 
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In order to find the renormalization constants, we need not the entire exact expression (6.13), rather we need its 
UV divergent part. The latter is proportional to a polynomial built of four facto rs dO at a single spacetime point x. 
The n eeded four gradients have already been factored out from the expression (6.11): two factors from the vertex 
( 3.12 ) and two factors from the ordinary vertices (2^). Therefore, we ca n neg lect the spacetime inhomogeneity of the 
gradients and replace them with their values at the point x. Expression (6.13) can therefore be written, up to an UV 
finite part, in the form 



where we have denoted 



dX3 / dX4 {di^6i^{x)ei {x3))o {di^0i^{x)9[^{xi))o {Vi^{x3,)Vi^{x4))o 



or, in the momentum-frequency representation, after the integration over the frequency. 



with Z3o from Eq. (1.3). Using the isotropy relations 



Pk/(A:). 



all possible permutations 



d(d + 2)---(d + 2n-2) 
the integral ( |6.16| ) can be reduced to the simple scalar integral 



m/(fc) 



Vk 



(6.14) 



(6.15) 



(6.16) 



(6.17) 



(6.18) 



with Cd from Eq. ( 2.11 ); the parameter m has arisen from the lower limit in the integral over k. The explicit answer 

for the quantity ( 6.16 ) is given in Appendix B. 

Contraction of the delta symbols with the first factors in Eq. (S.14) gives rise to various monomials built of four 

gradients of 9; up to total derivatives, they reduce to the operators from the family ( |6.7D . Then the function T^a\x; 9) 
from Eq. (6.9) in the one-loop approximation of the renormalized perturbation theory (i.e., to the first order in g) up 
to an UV finite part can be written in the form 



T^^){x-9)^F^ + ^-^{^ji/mY R^, 



(6.19) 



where /i has appeared from the relation Dq .9o^'o ~ 9^^^'^ ^-nd 



(6.20) 



are linear combinations of the monomials (6.7) with the coefficients A^p dependent only on d: 

(3rf2 + 22d + 36)Fi 2{d + i)F2 6F3 



Ri = 



i?2 =- 



i?3 



d{d + 2){d + A){d + Q) d{d + 2){d + A){d + Q) d(d + 2){d + 6) 

2{d + 5)Fi 4(rf + 3)F5 AFe 

'{d + 2){d + 6) d{d + 2){d + 4:){d + 6) {d + 2){d + 6)' 
12Fi (d^ + I0rf3 + 19^2 - Ud - 90)F2 



d{d + 2){d + 4:)(d + 6) 



d{d + 2){d + 4:)(d + 6) 



jd + 12)f3 
d{d + 2){d + 6) 



12Fa 



2{d^ + 8d^ + 10d-18)F5 



4.Fa 



d{d + 2){d + 6) d{d + 2){d + 4){d + 6) {d + 2){d + 6)' 



(d2 + 18d-t-48)Fi 



2{2d + 9)F2 



(d^ +4:d- 15)F3 



2d{d + 2){d + A){d + 6) d{d + 2){d + A){d + 6) 2d{d + 6) 

2Fi 2{d^ + 6d + 6)F5 2{d + A)F6 

^{d + 2){d + 6) ~ d(d-h2)((i-h4)((i-F6) (d 2)(d 6) ' 
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Ra 



2{d + 3)Fi {d + 3){d^ + 8d+U)F2 

d{d + 2){d + 4:){d + 6) ^ 2d(d+2)(d + 4)(d + 6) 
(d3 + - 14d - 36)^4 2((i + 3)(d+5)F5 



3Fs 



d(d + 2)(d + 6) 
2Fe 



i?5 = - 



i?fi — 



2d(d + 2)(d + 


6) 


d{d 


+ 2){d + A){d 


+ 6) 


(d + 2)(d + 6)' 






,{d' + 


lOd^ + 27d + 




3F3 


fi(d + 2)(d + 4)(<i 


+ 6) 


d{d 


-f 2)(d + 4)(d 


+ 6) 


d{d + 2){d + 6) 


2(d + 3)F4 


{dU 


- lOd^ + 


18^2 - 5M - 


114)^5 


12F6 



d(d + 2)(d + 6)' 



(d2 + i8d + 48)Fi 



(2d + 9)F2 



Ad{d + 2){d + A){d + 6) 
[d^ 



{d^ + 6d + 6)f3 

d{d+2){d + A){d + 6) ^ 2d(d + 2)(d + 6) 



ed + 6)F5 



(d + 2)(d + 6) d(d + 2)((i + 4)(fi + 6) 



+ 6rf2 - llrf-42)F6 
2d{d + 2){d + 6) 



(6.21) 



The constants Zap are found from the requhement that the functions (6.E) for the renormahzed analogs of operators 
(3.7), defined by the relation — ZapF^', be UV finite, i.e., be finite for e ^ 0. In the MS scheme this gives: 



Zcp - Sap + gCd Aap/e + 0{g^) (6.22) 

with the coefficients Aap from (6.21). For the matrix of anomalous dimensions (|6.1[ ) at the fixed point (4.9) one has 
7*^ = —g^CdAafj + O(e^), and for the matrix of critical dimensions from Table | and Eq. ( 4.11 ) one obtains 



= 2e6aP + lap- 



(6.23) 



Critical dimensions associated with the family of operators ( |6.7[ ) are given by the eigenvalues Aq of the matrix ( |6.23 ) 
In particular, for d = 3 one obtains: 



Ai w -0.55e, Aa « 0.68e, A3 « l.le, A4 ~ 2Ae, A5 = 8e/3, Ag = 3£. 
For general d, only one of the eigenvalues is found analytically: 

(d+l)2 



A, 



(d2-3) 



(6.24) 



(6.25) 



(we recall that all these eigenvalues have corrections of order O(e^) and higher). The critical dimensions Aq, as 
functions of d are presented in Fig. ^. They are always real, except for the pair A2,3 which becomes complex 
conjugate in the interval 4 < d < 5 (in Fig . H th e real part is shown). One can also see that for all d, exactly one 
of the dimensions, denoted by Ai in Eq. (6.24), is negative, and the others are positive. Existence of a negative 



dimension implies that the fourth-order structure function in model (1.1)-(1.3) exhibits the inertial-range anomalous 
scaling; for small mr it has the form 



^4(r) 



2e 



^Aa{mrY 



(6.26) 



see Eq. ([3.2|). The dots stand for the corrections of the form (mr)^+''('^^ and higher, which arise from the operators 
including more derivatives than fields, and possible anisotropic contributions, related to nonscalar operators. The 
leading term, singular for mr — > 0, is determined by the negative dimension Ai. 

The dimensions diverge for d \/3 as a result of the divergence of in Eq. (4.£). For d 00, they simplify 
and form three groups which tend to 0, e and 2e. More precisely, in the first group there are two dimensions with 
the behavior A = ±2^/2e/d + ©(l/d^). This means that for d = 00, the anomalous scaling of ^4 vanishes: the 
phenomenon known for the scalar Kraichnan model |^ and questioned for the NS turbulence In the next 

Section we shall see that the simplification of the exponents and vanishing of the anomalous scaling for d 00 also 
holds for the higher structure functions. 

Although the above expressions for the eigenvalues are well-defined for any d > a/3, low integer dimensions require 
special care because of additional linear relations between the operators. For d = 3, there are two such relations: 
^5 — $1/2 — $2 + 2<i>4 = and $6 — $3/2 = 0. For d = 2, one more relation arises: $4 — $i/2 — 0. (It is also 
noteworthy that the derivative on the right-hand side of Eq. ( |6.8| ) in two and three dimensions vanishes identically.) 

Using these relations, one can check that two basis operators for d = 3 and three basis operators for d = 2 
vanish. Therefore, the corresponding eigenvalues are in fact meaningless and should be discarded (in particular, this 
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happens with the dimension (6.25)). The eigenvalues that survive for d = 2 and 3 belong to the three and four 
lowest branches in Fig. ^, respectively (they are denoted by the thick dots). In particular, the dangerous operator 
remains nontrivial, so that our model exhibits the anomalous scaling also in two and three dimensions. We recall that 
similar behavior was demonstrated by the zero-mode solutions: all rf-dimensional expressions for the exponents have 
well-defined limits for d 2, but a part of them becom es in fact spurious owing to the vanishing of the corresponding 



amplitudes; see the discussion in the end of Sec. Ill B 



C. Scalar operators of the form {86) " and the anomalous scaling of S2n 

The leading terms of the inertial-range behavior of a higher-order structure function S2n are related to the scalar 
composite operators of the form {39)^^ with < A: < n. In the previous section we have established the anomalous 
scaling behavior of the fourth-order structure function, as a result of the existence of dangerous operator with k = 2 
in the corresponding OPE. Then the probabilistic inequalities allow one to show that all the higher-order structure 
functions are also anomalous, the leading term of the inertial-range behavior for the function S2n is given by an 
operator with k = n, the number of dangerous operators is necessarily infinite, and the spectrum of their dimensions 
is not restricted from below. 

Let A„ be the dimension of the operator that gives the leading contribution to the OPE for the function S'2„ , so 
that oc D(^"r"*^^~'^) (mr)'^" . It is well known in the probability theory that |(a;")|^/" is a nondecreasing function 

of n for any random variable x. Taking x — [6',.(i,x) — 6',.(t, x')]^ we find that S'2/1" oc Dg^^r'^^^^^ (mr)^"/" is a 
nondecreasing function of n, and so is the ratio |A„|/n [we recall that A„ is negative and (mr) is small]. This proves 
all the above statements. 

In principle, calculation of the critical dimensions related to the family (90)^" for any given n is a purely technical 



problem, and the formulas (6.11)-( 3.1(;) remain valid in the general case with obvious alterations. In practice, however, 
this problem appears very cumbersome, in particular, because the number of relevant operators increases rapidly with 
n. The situation simplifies for large d, and below we restrict ourselves with the zeroth and first terms of the 1/d 
expansion. To avoid possible misunderstandings, it should be stressed that we deal with the 1/d expansion of a critical 
dimension A in its 0(e) approximation, that is, the 1/d expansion of the coefRcient A'--^\d) in the representation 
A = eA(i)(d) + 0(e2). 

Despite this simplification, no explicit analytical result is available for general n. Below we only present the results 
for the critical dimensions of the families (96')^" with n < 6; the detailed derivation is given in Appendix B. 

The first two terms of the 1 /d expansion for such operators have the form A*^^) (d) — 2k + A^^^^ /d + 0(l/c?^), where 
A; = 0, 1, . . . , n and A'-^-'^' are numerical coefficients independent of the parameters e and d. It is clear that for large d, 



dangerous operators can only be present in the subset with fc = 0. As already mentioned in Sec. VI B, in the family 
with n = 2 there are two such operators with 

A(") = ±2^/2. (6.27a) 
In the family with rt = 3 there are three such operators with 

A'") = -9.674, -0.973, 7.647. (6.27b) 
In the family with n — A there are five such operators with 

A(") = -20.617, -7.783, -1.018, 3.883, 14.534. (6.27c) 
In the family with n ~ 5 there are seven such operators with 

A(") = -35.589, -18.660, -8.700, -2.960, 2.780, 10.674, 23.455. (6.27d) 
In the family with n — 6 there are eleven such operators with 

A(") = -54.572, -33.612, -19.554, -13.834, -12.815, -4.908, 

3.839, 4.828, 9.681, 19.552, 34.395. (6.27e) 

These results confirm and illustrate the general picture outlined above: in the set of operators with k < n, the 
most dangerous operator (that is, the operator with the lowest negative dimension) belongs to the subset with k = n, 
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and its dimension A„ < decreases faster than linearly with n. The results (3.27) are illustrated by Fig. ^ (only 
the negative dimensions are shown). It suggests that the dimensions form a set of monotonous branches, denoted by 
dashed lines. The solid curve corresponds to the well-known expression for the scalar Kraichnan model |^,|| in the 
same 0{e/d) approximation: A„ — —2n{n— l)e/d. For all n, it lies below the lowest vector branch (the scaling in 
the scalar model appears "more anomalous" ) ; the deviation between the scalar and vector cases becomes stronger as 
n increases, although the ratio of the dimensions approaches unity. 



D. Tensor operators and the scaling of S2 in anisotropic sectors of arbitrarily high orders 



In this Section, we apply the RG and OPE approach to the higher anisotropic sectors of the model (l.l)"( |l.3[ ). 
We shall concentrate on the second-order structure function, for which nonperturbative results c an, in principle, be 
derived for arbitrarily high values of the parameter I from the exact Dyson- Wyld equation (see Sec. [II B ). This allows 
one to identify the solutions of the zero- mode equations, discussed in Sec. [II, with definite composite operators and, 



in principle, to calculate the corresponding amplitude factors. Using the OPE technique, we derive explicit analytical 
expressions for the leading exponents in all anisotropic sectors to the order 0(e) in d dimensions. Furth erm or e, w e 
present additional nontrivial exponents which do not appear in the inertial-range behavior of the model (l.l)-( pT3| ), 
but will be activated (and can determine leading terms in anisotropic sectors!) if the anisotropy is introduced by the 
velocity field (like in Ref. ) and not only by the large-scale forcing. 

The analysis of the anisotropic sectors for higher-order functions using the OPE is extremely cumbersome but, in 
a sense, purely technical problem; we shall briefiy discuss it in Sec. VI E . 

According to the general rules (see Sec. 0), the leading terms in the sector I — 2 are determined by the second-rank 
operators built of two gradients; up to derivatives, there are two such operators: 



= inr \^^0k^Jek\ , F2 - inr \^ke^^ke^ 



(6.28) 



where JTZV denotes the irreducible part; cf. Eq. ( 3.11 ) in Sec. [II B (here and below, we use Fa and to denote 
different operators and their dimensions). 

We omit the one-loop calculation of the corresponding 2x2 matrix of critical dimensions, which is similar to the 
calculation discussed in Sec. VI B for scalar operators, and present only its eigenvalues, i.e., the critical dimensions 
of operators ( 6.28| ): 



Ai,: 

I n rep resentation 



d^ 



2d-8± y/{d + 4)(d5 + 2d4 - 7^3 _ 4^2 + 8d + 16) 



2(d2-3)(d + 4) 



+ 0(e2). 



(6.29) 



4), they give rise to exponents 2 — e -f Ai^2, which agree with the special case Z = 2 of expressions 



), (3.23), obtained in Sec. IIIB on the basis of the Dyson-Wyld equation. 



In two dimensions, the transverse vector field can be represented in the form 



(6.30) 



where e^^ is the antisymmetric Levi-Civita pseudotensor and tpi^) is some scalar function (stream function). Using 
the well known identity Cik^i s = S gSks — ^is^kj one can easily check that for c? = 2, the operators (3.28) coincide. 
Only one of the dimensions ( 6.29 ), namely, A2 — e, corresponds to a nontrivial basis operator (6.2) and remains 
meaningful, while the other, Ai — 3s, corresponds to a vanishing basis operator and should be discarded; cf. Fig. |^ 
for d = ; = 2. 

The leading terms in the sector I = 4 are determined by the fourth-rank operators, obtained from the monomial 
Fijkn = mfP \di9jdk0ri\ by all possible permutations of its tensor indices. It turns out that there are only three 
different critical dimensions, associated with these operators. The corresponding basis operators (6.2) possess different 
symmetries and therefore can be written down without calculation of diagrams. One of them is the fully symmetrized 
operator. 



Fs — Fijkn + Fjink + Fikjn + Fk 



jikn 



F, 



iknj ■ 



The others can be constructed as follows. The monomials can be split into three groups of two operators each, 

-^1 Fikjn 1 F'2 Fkinj ; 

symmetric with respect to the simultaneous exchange of the indices within the pairs {ij} and {kn}, 
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-^3 F'ljkn-^ -^4 ^jink: 

symmetric with respect to the exchange of the indices within the pairs {ik} and {jn}, 

symmetric with respect to the exchange within the pairs {in} and {jk}. 

The operators can mix in renormahzation only within the groups with the same symmetry, so that the set Fi, F2, 
Fs is closed wit h re spect to the renormahzation, and so are the sets F3, F4, Fs and F5, Fq, Fs- In the first set, the 



basis operators (6^) are Fs (it is fully symmetric and no other operators can admix to it), Fi — F2 (it is antisymmetric 
with respect to the exchange of the indices within the pairs {ik} and {nj}), and Fs — 3Fi — 3F2. The latter operator 
is symmetric with respect to the exchange of the indices within the pairs {ik} and {nj}, and for this reason it cannot 
mix with the second basis operator; it is not fully symmetric and cannot admix to Fs- It remains to note that the 
contraction of the operators Fi — F2 and Fs — 3i^i — 3F2 with any constant vector with respect to all four indices gives 
zero, while analogous contraction of Fs remains nontrivial. This means that Fs cannot admix to those ope rators in 
renormahzation, and the above construction indeed gives three independent basis elements of the type (|6.2[). 

The explicit one-loop calculation confirms this conclusion and gives three different critical dimensions, corresponding 
to the basis operators Fs, Fs — 3Fi — 3F2 and Fi — F2, respectively: 

(d + 2)((f +4:d-9) , 



(6.31b) 



(d2-3)(d + 4) 



(6.31c) 



with corrections of order 0(£^) and higher. 

The remaining independent basis operators can be obtained by permutations of the indices and can be chosen in 
the form F3 — F4, F^ — Fq and Fs — 3F3 — 3F4. At first glance, it seems that one can add another independent element, 
Fs — 3F5 — 3i^6i but in fact it is equal to the su m of the operators Fs — 3Fi — 3F2 and Fs — 3^3 — 3F4 up to the 
minus sign. Therefore the dimension Ai in Eq. ( 6.31 ) is unique, A2 has two- fold degeneracy and A3 has three- fold 
degeneracy. 

Although all three dimensions in Eq. (3.31) make sense, only Ai appears on the right-hand side of expansion 
(5.2). Indeed, both the diagrammatic analysis and dimensional considerations show that the coefficients in Eq. (^.l|), 
corresponding to operators Fij^n, do not involve the function {69)o from Eq. ( 2.4b| ). Therefore, they do not depend 



on the vector n and their tensor indices are carried by the Kronecker delta symbols or vectors r. The contraction of 
any irreducible operator with the delta symbols always gives zero; the contraction with the components of a single 
vector r in expansion ( ^.1[ ) "kills" the basis ope rators of the form Fs — 3Fi — 3F2 and Fi — F2 (see above) . 

Thus the only contribution to expansion ( |5.2| ) comes from the operator Fs, and the expression 2 — e -t- Ai should 
be identified with the leading exponent 74 from Eq. ( 3.21 ) in Sec. Ill; they indeed agree in the order 0{e). 

Th e remaining dimensions A2,3 will be activated (and can determine leading terms on the right-hand side of Eq. 
(5^) in the I = 4 sector) if the vector n appears in the corresponding coefficients Cf in expansion (^!1|). This can 
happen if the anisotropy is introduced by the velocity field (like in Ref. p^) and not only by the large-scale forcing. 
However, it should be noted that in such a case the dimensions become nonuniversal through their dependence on 
the anisotropy parameters, and the expressions like (6.31) give only the zeroth order (isotropic) approximations; cf. 
Ref. for the scalar case. 

Now let us turn to general I. The relevant operators are the Ith rank tensors built of two fields 9 and minimal 
possible number of derivatives: 



Up to total derivatives, m ono mials (6.32a) are invariant with respect to the shift 6 
can appear in expansion (| 



.32a) 



const, so that their dimensions 



.2[ ). In the above form, all symmetries of operators (3.32a) are obvious: they are symmetric 
with respect to the permutation of the indices {ii,ii} an d any permutations within the subset {«2 • • -ii-i}. Thus 
the total number of different monomials, obtained from (6.32a) by permutations of the indices, equals l{l — l)/2. 
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However, there are only three different dimensions related to them. Indeed, the co unterterm to the monomial (6.32a 
necessarily has the same symmetries and therefore can include, along with ( |6.32a ) itself, two more structures: 



SymlTZV [di^di^di^ ■ ■■di^_-^l 



SyminV [Oi^di.di^di^ ■ ■■di,_-^9i,, 



(6.32b) 



where Sym denotes the symmetrization with respect to the perm utation of the pair {ii,ii} and any permutations 
of the indices {12 • • -ii-i}. Th us t he set of three operators (6.32) is closed with respect to the renormalization, the 
corresponding basis operators (6^) are their linear combinations and determine three different dimensions Ai^2.3; all 
the other basis operators are obtained by permutations of the indices and give rise only to the same dimensions. 

For general I, even the one-loop calculation is rather difficult because individual contributions in the counterterms 
to the polynomials (3.32) co ntain powerlike UV divergencies in addition to logarithmic ones. In contrast with the 
calculation discussed in Sec. VI B, one cannot neglect the spacetime inhomogeneity of the fields 9 in the diagram; 
in other words, one cannot neglect the dependence of its integrand on the external momenta p and should expand 
the integrand up to the terms of order This expansion gives rise to the terms that diverge for A ^ 00 as 

some positive powers of the UV cut-o ff A. However, all such terms contain delta symbols and cancel out when all 
contributions in irreducible structures (S.32) are taken into account; the result is fini te at A 00, contains a first-order 
pole in £, and reduces to a linear combination of the three structures ( 6.32a ) and ( 6.321 ). 

We omit the details of this cumbersome calculation and give only the result: 



A„ = (Z - 4) + 



S Xn 



{<P - 3)(d + 2l~ 6)(d + 21- 4)(d + 21-2) 



(6.33) 



where 



xi^d^ + MH + liSf + 10(fP + dl"^ - I2d^ - 5'Sd^l - GOd^l^ - 6dl^ 
+92d^l - 29dP - 44/^ - 24^^ + 158d/ + 214^^ - 156d - 



+ 21^- 



47d^ 
192, 



X2 = id + 2l- 4)(d + 21- 2){d^ + 2d^l + df - Qd^ - 5dl + 2f + 'id- l&l + 30), 



xj. = {d + 2l- 2){d'^ + UH + + dp - lOd^ - 2bdn - 5df + 2P + 27d^ - 6dl - 2QP + 30d + 92^ - 96). 



For I — 4, the results ( |6.31| ) are recovered. Like in the case I = 4, the basis operator ( |6.2D that possesses the dimension 
Ai is symmetric with respect to all possible permutations of the full set {ii • • -i;}; only this o per ator survives the 
contraction with the coefficients Cp in the operator product expansion (5T) for the model (p~l|)-(1.3). Therefore only 
Ai appears on the right-hand side of Eq. ( |5.lD and determines the leading exponent for the Ith anisotropic sector: 
7; = 2 — e + Ai. For all I and d, this recovers the result ( 3.21 ) obtained in Sec. [H on the basis of the Dyson- Wyld 
equations. 

Like in the case I — 4, the remaining dimensions A2,3 in Eq. ( 6.33| ) are activated and appear on the right-hand side 
of Eq. (5.2) when the anisotropy is introduced by the velocity field. 

We recall that the general form of the exponent for a given ^ > 2 is 7/ = ^ — 2 -I- 2fc -I- 0(e)] see Sec. [II B . From the 
OPE viewpoints. A: = corresponds to the I th rank irreducible operator built of two fields 9 and I — 2 derivatives, 
symmetric in all indices (see above). The operators with A: = 1, 2, . . . can be obtained in two ways: one can add k 
Laplacians to the operator described above, or one can add (fc — 1) Laplacians, two derivatives with free indices, and 
contract the indices of the 9 fields. Therefore, for general d the leading exponent 7/ — I — 2 + 0{e) is unique, while 
for all /c = 1, 2, . . . there are two correction exponents o f the form 7; = Z — 2 -I- 2A; + 0(e). In two dimensions, these 
two possibilities coincide, see the discussion below Eq. (3.2£), and only one exponent exists for any fc = 0, 1, 2, . . .. 

For I — 0, the general form of the exponent is 7 = 2fc + 0(e); the fir st two solu tions , 7 = and 7 = 2 — e, arc known 



exactly both from the RG and the Dyson- Wyld equation; see Sees. HI A and VIA. For any k and rf, the solution 
is unique: for the scalar operators, the indices of the fields are contracted, and the operator with 7 = 2fc + 0(e) 
necessarily reduces to the unique form 9iA''9i. 

This picture is in a full agreement with the results obtained in Sec. Ill for general d and / on the basis of the exact 
Dyson-Wyld equation for the pair correlation function; see also Rcf. [26[ for d = 3 and j < 10. 



E. Higher-order structure functions in the higher-order anisotropic sectors 

Let us briefly discuss the scaling behavior i n the anisotropic sectors for the higher-order structure functions. The 
RG and OPE analysis given in Sees. and |VID| can directly be extended to the general case. It shows that the 
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leading exponents 7„i in the I th sector of the 2n-order structure function, S2n oc Pi (z) r'^"-' , are determined by the / th 
rank tensor operators with k < 2n fields 9 and minimal possible number of derivatives; the operato rs w hich contain 
the field 9 without derivative, or reduce to total derivatives, give no contribution to the expansion (5.2) and should 
be discarded. 

The practical calculation of the critical dimensions of the operators with large I or n is a difficult task, as one could 
see already on the example with I = and n = 2. However, in the zcroth order of the e expansion, some important 
information can be obtained without calculation, just by the analysis of the form of relevant operators. 

One can easily see that for I < An, the leading exponents are determined by the I th rank tensor operators of the 
form {89)^'^, with I free and An — I contracted indices. In the 0(1) approximation, the exponents themselves are equal 
to the number of derivatives entering into the operators: 7„/ = 2n + 0{e). 

For / > An, the relevant operators necessarily contain more derivatives than fields, and for the leading exponents 
one obtains: 7„; = I — 2n + 0{e). 

One can thus conclude that for higher n, the general picture remains the same as for S'2: each anisotropic sector 
possesses its own set of scaling exponents; the leading exponents obey hierarchy relations at least for small e and 
I > An; they grow with I without bound. Of course, there can be several exponents for given I and n; we recall that 
there are six exponents of the form 720 = 4 + 0(e) for I = and n = 2 (see Sec. VI C). In order to identify the unique 
leading exponent within a family with the same zeroth-order value, or to verify the hierarchy relations for I < An, 
one should perform the 0{e) calculation for the relevant families of operators. This cumbersome task lies beyond the 
scope of the present paper and will be discussed elsewhere. 



VII. CONCLUSION 



We have studied the inertial-range scaling behavior in a model of the passive vector quantity advected by a self- 
similar white-in-time Gaussian velocity field, with the large-scale anisotropy introduced by a random forcing. In 
two respects, the model is closer to the real Navier-Stokes turbulence than the famous scalar rapid-change model: 
nonlocality of the dynamics and mixing of the composite operators that determine anomalous exponents. 

The incompressibility condition for the advected field and the pressure term in the diffusion-convection equation 
make the dynamics nonlocal. This raises the question of rcalizability of the zero-mode solutions, that is, convergence 
of the integrals in the equations for the correlation functions on power like solutions, and consistency of nonlocality 



and the existence of infinite families of scaling exponents [g6 36 1 . The detailed analysis of the exact integral equation 
satisfied by the pair correlation function has shown that the general picture of the inertial-range scaling is essentially 
the same as in the scalar and magnetic variants of the rapid-change model. Namely, each anisotropic 

sector is described by an infinite set of scaling exponents, with the spectrum unbounded from above. The leading 
exponents in each sector are organized in the hierarchical order according to their degree of anisotropy, with the main 
contribution coming from the isotropic sector in agreement with the hypothesis on the restored local isotropy of the 
fully developed turbulence The leading exponents themselves grow without bound with the degree of anisotropy, 
in disagreement with the idea of the window of locality . 

The integral operator entering into the equation for the pair correlation function in the momentum space converges 
on the powerlike solutions with the leading exponents in the I = and I = 2 sectors, but formally diverges on 
powerlike solutions with subleading exponents and leading exponents in the sectors with I > A. However, correct 
analysis of convergence here requires the knowledge of the behavior of the full solution beyond the inertial range, 
where it no longer reduces to a sum of power terms. It turns out, that natural assumptions about the form of the 
solution allow one to perform certain subtractions in the integrals which make them convergent. One can use the 
formal rules of analytical regularization and simultaneously omit the subtracted terms to obtain correct answers for 
the convergent integrals with proper subtractions. Moreover, the realizability of these solutions is also guaranteed by 
the RG approach, where they are identified with the contributions of certain composite operators in the corresponding 
operator product expansions. Therefore, such exponents indeed appear in the full solution in the inertial rang e. 

These conclusions are in agreement with the recent analysis performed in Ref. for the model (!.!)-( 1.3) in the 
three-dimensional coordinate space, although the analysis in the momentum space appears rather different; see also 
Ref. for the general vector model. Furthermore, the RG and OPE techniques confirm this picture and extend it 
to the higher-order correlation functions. 

The second aim of the paper has been the analysis of the anomalous scaling of the higher-order even structure 
functions S2n- Owing to the conservation of the "energy" 9'^{x), the second-order function appears nonanomalous 
with the simple dimensional exponent: S'2 oc r^~^. The anomalous scaling reveals itself on the level of the fourth-order 
structure function. In contrast with the scalar case, where the leading anomalous exponents were identified with the 
critical dimensions of individual composite operators in the corresponding OPE pf , the vector nature of the advected 
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field in our model leads to mixing of operators. In particular, the inertial-range behavior of the function S4 in d 
dimensions is given by a set of six close exponents, determined by eigenvalues of the matrix of critical dimensions 
for a set of six operators. One of the dimensions is negative ("dangerous operator") and gives rise to anomalous 
scaling. The number of relevant operators increases rapidly with the order of the function; they have been calculated 
in a controlled approximation (small e and large d) for the higher-order functions up to 5*12. The latter involves as 
many as sixteen negative exponents, ten of them coming from the lower-order functions, and a multitude of positive 
exponents which are small and therefore close to the negative ones for small e. The probabilistic inequalities prove 
that all the higher-order structure functions are also anomalous, and the total number of dangerous operators in our 
model is infinite, with the spectrum of dimensions unbounded from below. 

Since the mixing of operators is a manifestation of the vector nature of the advected field, there are a few general 
conclusions regarding the real NS turbulence, which one can draw from the analysis of model (1.1)-(1.3). 

It was demonstrated recently that a careful disentangling of contributions from different anisotropic sectors is 
ubiquitous in the analysis of experimental data on the real turbulence, because it allows one to properly identify 
scaling exponents in the situations, where the standard treatment reveals no s caling behavior at all; see the discussion 



Rcfs. ||32|-|35| and references therein. The example of the model (IT)-(l^) shows that even in the isotropic sector 



or for the ideal isotropic turbulence, correlation functions are represented by infinite sums of powerlike terms, and the 
number of close terms grows rapidly with the order of the correlation function. Although these corrections die out 
in the formal limit L — > 00, and a pure powerlike behavior with the leading exponent sets in, in practice it may be 
obscured by such corrections: the subleading exponents can be very close to the leading ones and much more important 
than the leading terms from the higher anisotropic sectors. This might result in imaginary nonuniversality of the 
inertial-range exponents or deviations from a pure scaling behavior, which increase with the order of the correlation 
function. Therefore, reliable analysis of the inertial-range scaling necessarily requires some theory for the correction 
exponents. 

In theoretical models, anomalous scaling is usually explained by the so-called intermittency phenomenon. Within 
the framework of numerous models, the anomalous exponents are related to the statistics of the local dissipation 
rate or to the dimensionality of fractal structures formed by small-scale vortices in the dissipative range; the detailed 
review and bibliography can be found in Rcf. [0. As a rule, those theories predict simple analytic formulas for the 
dependence of the anomalous exponents on n, the order of the structure function. Although such formulas can provide 
a very good fit for the experimental results, the experience on the rapid-change models suggests that they cannot be 
absolutely correct. Even for the scalar model, the n-dependence of the anomalous exponents changes as the order of 
the £ increases j^-p^. Moreover, for our vector problem the n-dependence of the exponents can hardly be given by 
a single explicit formula (except probably for the case d = 2, where the model can be mapped onto a nonlocal scalar 
problem): the relevant families of operators are completely different for different n, so that each function Sn requires 
special analysis. 

Recently, a systematic perturbation theory for the anomalous exponents in the NS turbulence was proposed, where 
the role of a formal expansion parameter is played by the anomalous exponent for S2, assumed to be small but 
nontrivial jisj. In contrast with the magnetic variant of the rapid-change model [^7y28[ and the general "^-model" 
with a stretching term p4| , the second-order structure function in our vector problem ( |l.l| )-(1.3) is nonanomalous, and 
the perturbation theory of Ref. |Q would be impossible here. Like in the scalar model |§-5U|, this is a consequence 
of the energy conservati on or , in the field-theoretic language, of vanishing of the critical dimension of the local 
dissipation rate; see Sec. VI A for the vector and Ref. |^ for the scalar cases. The vanishing of the critical dimension 
of the dissipation rate at the physical value of e is also characteristic of the NS case , which raises serious doubts 

about the existence of the second-order anomaly and the possibility of the corresponding perturbation theory in the 
NS turbulence. 

The analysis of the inertial-range behavior essentially simplifies as 0? — s- cx). Our model has no finite "upper critical 
dimension," above which anomalous scaling vanishes (see Ref. for a recent discussion of that concept). Like in 
the scalar case and, probably, in the NS turbulence the anomalous scaling disappears at d — 00, but it 

reveals itself already in the 0{l/d) approximation. The anomalous exponents can be calculated within the double 
expansion in e and 1 /d. Along with the results for the scalar rapid-change model, this confirms the importance of 
the large-(i expansion for the issue of anomalous scaling in fully dev elop e d tu rbulence. 

Although our analysis has been confined with the linear problem (^)-(|l^), which has only restricted resemblance 
with the real fluid turbulence, some of the results can be extended to the case of the vector passive field advected by 
the NS field, or the nonlinear NS equation itself with some classes of random forcing. These questions lie beyond the 
scope of the present paper. Detailed exposition of the RG approach to the NS problem and the bibliography can be 



found in Refs. |16 V^; the renormalization of composite operators and the concept of the operator product expansion 
are also discussed in Refs. | |l5| , ^ , ^ , ^ , p7||4S|j50t . In particular, critical dimensions of tensor composite operators in 
the stirred NS problem were calculated in Rcf. (see also Sec. 2.3 of [|l^); they demonstrate the same hierarchy 
as their counterparts from Sec. VI D. 
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We believe that the framework of the renormalization group and operator product expansion, the concept of 
dangerous composite operators, exact functional equations, and the e and 1/d expansions will become the necessary 
elements of the appearing theory of the anomalous scaling in fully developed turbulence. 



ACKNOWLEDGMENTS 



The authors are thankful to M. Hnatich, J. Honkonen, A. Kupiainen, A. Mazzino, P. Muratore Ginanneschi, 
S. V. Novikov and A. N. Vasil'ev for discussions. The work was supported by the Russian Foundation for Fundamental 
Research (Grant No. 99-02-16783) and the Grant Center for Natural Sciences of the Russian State Committee for 
Higher Education (Grant No. EOO-3-24). 



APPENDIX A: EQUATIONS FOR THE EXPONENTS IN THE HIGHER ANISOTROPIC SECTORS OF 

THE PAIR CORRELATION FUNCTION 



It was shown in Sec. Ill B| that the equation for the exponents 7; in I th anisotropic sector can be written as 
detjC^'jil = 0. For I ~ 2, the matrix elements C^'^ were given in Eq. ( 3.19| ). In general, the matrix C^'^ is 
symmetric and its elements are finite linear combinations of the integrals Ii = Ii -\- 2 J with J from ( 2.13 ) and 

In = Iniji — 2,£ -I- 2) from ( |3.3| ), with n as high as 1/2 + 2. In this Appendix, we present the coefficients C^'j, for 
higher values of Z up to ^ = 12. Then equations for 7/ can be written in a straightforward way. Below we denote 

dk^s = {d + k){d + k + 2)...{d + s-2)(d + s) and dk = (d + k). 

= ^^2,4/4 ~id + 2){d^ + 5d- 2)h + (d' - l)(2d + 5)/2 -{d- lf{d + l)Ji, 

rr(4) 



-^2,4/4 + 2^2,4/3 - {d+l){2d+b)h + {d^ - l)h, 

d2Ah - \d2Ad + '^)h + -^^7^ (rf' + lOd' + 24d+ 12)/2 - ^{d^ - l)(d + 4)7i, 



12 



12 



'-'11 



(d2 + iod+l)^ d3(d' + 6d-5)^ di3(d-l)(3(i+13)^ di3(d-l)2~ 

-h - 3 ; 73 -I ; 12 '■ ; il. 



rfs 



di.i 



di.i 



(6) 



^ (3d +19)^ , d3(4d + 25)^ di,3(3d+13)^ ^ d^^d-l) ^ 
; -'4 H ; -13 ; ^2 H ; -ii, 



^6,8 



di.s 



^4.8 



^ ^ (^2 + l4d + 43)^ ^ d3(d^ + 30^2 + 216rf-h 430)^ ^ ^1,3(^2 + 15^2 + gg^ _^ 85)^ 

^ 3(i8 30c?6,8 ^ 15c?4,8 ^ 

diAd-l){d^ +M+10)- 



30rf4,i 



.8) r (c«' + 15rf + 8)^ d^[2d-l){d+ll) d3,5(3d2 + 23d-22) di,5(d - l)(4d + 25) 

L-il = ^6 -\ + ^ -\ J4 — ^ -\ -'3 H -\ ^2 - 

"12 "10,12 "8,12 "6,12 



dl,5(d- 1)2 



-'6,12 



C- 



(8) 
12 



(2d+17) d5(7d + 64) d3^(7cr 
-H + ^ ; ; -'4 H ; — 



du 



)^ di,5(4d + 25)^ , di.5(d-l)- 
-J3 ; -'2 H ; J^li 



de,i2 



^6,12 



a 



(8) 
22 



^ (d^ + 24d-H16) ^ ^ rfs (^^ + 58^2 + 748d -f 2632) ^ ^3,5 (3d^ 90^^ -f 788rf + 2072) ^ 
4di2 ^ 56(iio,i2 56^8,12 ^ 



di 5(3^2 + 50d -I- 196) ^ di 5(d - l)(d2 + 12d -f- 28) ~ 
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22 
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15dio,i6 ^ 45^8,16 ^ 904,16 ^' 

(d2 + 25d + 34)^ „d9(2d2 + 39d_ 1) ^7 9(3^2 + 49d - 28) ^ d^ g{A(f + bbd - A7) ^ 
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'12,20 
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132dio,20 132dio,20 
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APPENDIX B: CALCULATION OF THE CRITICAL DIMENSIONS OF THE OPERATORS {def" FOR 

LARGE d 



T he cr itical dimensions of the scalar composite operators of the form (90) 2" were presented and disc ussed in 
Sec. VI C in the first nontrivial orders, 0(e) and 0(e/d), of the double expansion in e and 1/d; see Eqs. ( [6.271) . Below 
we give the derivation of those results. 

The UV divergent part of the diagram from Eq. ( |6.9| ) with a composite operator {dd)'^" is given by the formulas 
(S.9)-(6.2C) for n — 2; see Sec. VIB. With obvious alterations, they apply to the case of general n and can be 
summarized as [the 0{g) approximation in the renormalized variables]: 



2e 



where we have used the relations ( |l.4D and ( |6.18D and denoted 



(Bl) 



(B2) 



The integration over the unit sphere in d-dimensional space in Eq. (B2) can be explicitly performed using the isotropy 
relations (6.17): 



T, 



ni2J3U,Jlj2j3j4 



+ 



^ <5'ii'i3 <^i2j2 '^i4j4 ^jljs ^ (^^i2j2 ^iiji ^ilisjljs + ^jlj3 ^1232 ^iiisiiji + '^i4j4 '^Jlj3 '^ni2i3i2^ ~^ 

] ( ^jlj3 ^ili2i3Ui2ii + '^i2i2 ^ilhiijljiji + '^i4i4 '^il»2»3ili2j3 ) ^ '~i ^il 
'0,4 ^ ^ "0,6 



«2«3*4jlj2j3j4 



(B3) 



where do,fc = d(d + 2) . . . (d + fc); cf. Appendix A. 

Although the coefficients of the tensor (B3) behave as 0(l/d) for d — ^ 00, contractions of their indices ik can 
compensate the smallness (the indices jk are contracted with the factor dOdO in Eq. (Bl) and cannot be contracted 
with each other). Such contractions occur when the factor d^F{a) / dada contains the terms of the form Si-^^i^ Si^i^ or 
SiiU ^i2i3 (the third term, Si^i^ Si^i^, is forbidden by the transversality of the field 9). In particular. 
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Here and below, we use the notation = dj^9i^ . . . dj^9i^: cf. Eq. (|67|) in Sec. |VIB 



However, the both such contractions can appear simultaneously only in the term with do, 4 in Eq- (BS) and giv e 
an 0{l/d) contribution. The leading 0(1) terms arise from the contraction of the pair iii^ in the first term of ( |B3| ). 
The needed term of the form Si-^i,^ in the factor d^F{a)/dada appears in every differentiation of the block of F that 
contains the contraction of the indices of derivatives: 



(B5) 



Substituting Eq. 



into Eq. (Bl) gives the contribution 
gCd if^/m) 



2e 



gCd (p/m) 
2e 



be 



that is, the block is reproduced, and the counterterm to the operator Fa in the order 0(1) is proportional to 
the same monomial Fa. The number of such contributions equals to the number of the contractions between the 
derivatives in the monomial Fa, and we obtain: 



Aal3 = -^naSap + 0{l/d), 



(B6a) 



Taii 



-UaSaf) + 0{l/d), 



(B6b) 



Ai°) = (n-n,)e + 0(l/d). 



(B6c) 



where < nc, < n for operators from the family (90)^". The minimal possible value, Aq = 0, is reache d fo r Ua = n, 
that is, for the operators where all the derivatives are contracted only with each other. In the family (6.7), one has 
Ua = n = 2 for F2 and F5; = 1 (A^ e) for F3, F4 and Fg; = (A^ = 2e) for Fi. 
Let us turn to the calculation of the 0{\/d) correction to the results (B6). We write 



Ai°) + Ca/d 



with aL°^ fr om (|B6cD and numerical coefficients Cq determined by the relation 



Aet\A- -Ca/d\ = 0, 



(B7) 



(B8) 



where A is the matrix (4.11) for the family (96*)^" with a given n. 

Since in the 0(1) approximation the matrix A is diagonal with the diagonal elements Ai"'', all the nondiagonal 
elements of the matrix in Eq. (p8|) are of order 0{l/d), as well as its diagonal elements that correspond to the 



(degenerate) eigenvalue Aq°\ The diagonal elements that correspond to the eigenvalues different from A^'' are of 



(0) 



order 0{\). It then follows that the determinant in Eq. ( [B8| ) is of order 0(l/d ), where N is the degeneracy of the 
eigenvalue Aq^-* , and the vanishing of the full determinant in the leading approximation is equivalent to the vanishing 
of its X subdeterminant that corresponds to Ai"-* . 

This means that in the 0(l/c?) approximation, the equations for the coefficients Ca corresponding to different values 
of Aq°^ or, equivalentl y, ng , are independent, and these coefficients can be sought separately. 



with A^i*^ — 0, and below we confine ourselves with this family. For such operators, Ua = n (tensor indices of the 
derivatives are contracted only with each other, and the same holds for the indices of the fields), and they always can 
be represented in the form 



It is clear from Eq. (B6c) that for small 1/d, dangerous operators (Aq, < 0) can be present only among the operators 



F=(0l)"H02)^ 



(B9) 



where X]fe=i ^^fe ~ ^^'^ 4'k is the scalar operator that contains 2k fields 9 and cannot be represented as a product 



This operator necessarily reduces to the form [cf. Eq. (B4) for k = 1] 
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[SfcSl, S1S2, S2S3, ■ 



(BIO) 



Now let us collect all possible contributions of order 0{l/d). 



(i) As already mentioned above, an 0{l/d) contribution appears from the term with do, 4 in Eq- (B3), when the 
both derivatives in the vertex (6.12) act on 0i. In (Bl), this gives the contribution of the form 



2ed \mJ ■'^ 2ed \m 



that is, the operator F in Eq. (B9) reproduces itself, and the number of such terms is equal to ni. Therefore, the 
corresponding contributions to the matrices A^p and 7*^ in Eq. (B6) are diagonal and have the forms 



(Bll) 



ii) An 0{l/d) contribution appears from the term with do. 2 in Eq. (B3), when the both derivatives in the vertex 
( 3.12 ) act on any one of the n factors <i>^^ in F and produce the delta symbol Si-^i^. Each of these differentiations 



gives into Eq. (Bl) the contribution 



.||(ii)'x3.,«.a,..„ 



and the corresponding contributions into the functions (B6) are again diagonal: 

3n ^ 



2d 



s: * 3n 



(iii) The contributions into Eq. (Bl) from the first term of Eq. (B3) 

gCd(^iY d^F{a 



(-) 



The operation dj-^Oi^ djdai^i^ breaks the chain of contractions in 0^: 

dj^di^ dcj)k/dai^i^ = 2k [jiSi, siS2,S2S3, ■ 



, Sfc-l«l 



(B12) 



(B13) 



(B14) 



in the notation of (BIO). In the follow ing, the field with the index ji is not differentiated, since it does not belong to 
the operator F[a) in the vertex ( 6.12| ). 



The operation dj-^9i^ d/dat-^i^ acts either onto the factor (B14) or onto some other factor 0s in the operator (p3S|). 
In the latter case, another broken chain of the form (B14) appears; along with the first broken chain it gives rise to 
the unbroken chain with k + s elements, that is, (j)k+s- Therefore, this process gives rise to the counterterm 



f^(^)\4ks^,,s 
2ed \mJ 



(B15) 



for any pair of factors (ps in the operator (B9), and they determine nondiagonal contributions to the matrices 



The action of the operation dj^Oi^d/dai^i^ onto the chain ( B14 ) produces a number of different terms. One 
possibility is the breakdown of the chain of contractions of two kinds: 



[j 



isi, S1S2, 



, Sp-lJl,«lSp+l, . . . , Sk-lll\ 



(B16a) 



[jisi, S1S2, • ■ ■ , Sp_iii, jiSp+i 
The first variant is obviously (j)p(j)k-p, that is, the "decay" 

fc-i 



, Sfc-l«l • 



(B16b) 



(B17) 



The second variant gives (fc — 1) factors 4>k, and they give the diagonal contribution to the matrices (B6) 
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I 1 

SaAafj = 2^ X! "^kkik - 1) Sa/s 



(B18) 



fc=2 



with n/j and q from Eq. ( ]B9| ). 

Finally, the differentiation of the rightmost factor in Eq. ( B14| ) gives d^fc. This is the leading 0(1) contribution, 
but it also gives the 0{l/d) term after the substitution 



g^g, = — [l + -]+ 0{l/d 



The contribution to the matrices (B6) is also diagonal: 



s: A ^ s: <: * 



Collecting the contributions (§ll|), ( |B12D , (|B18|) and ( |B19| ) gives 



SAa/i = ^ (^ni - " + ^ nkk{k - ifj Sap, 



k=2 



(B19) 



(B20a) 



(B20b) 



fc=2 



for the total diagonal 0{l/d) contribution to the matrices (B6), while their nontrivial nondiagonal elements are 
determine d b y Eqs. ( B15| ) and (B17), with the summation in the former over all pairs of the factors (j)k4's in the 
operator (B9). 



Consider a few examples that illustrate the above algorithm and lead to the results announced in Sec. VI C 
For n = 2, the family of operators with n 



n = 2 consists of two elements, 
F^{(f>l 02}, 



and the equation (|BS| ) for the coefhcients Cq in the representation (B7) has the form (here and below we omit the 
subscript a and change the signs of the matrix elements so that they all become positive) : 



C 4 
2 C 



0. 



(B21) 



The solutions are C = ±2-\/2, as announced in Eq. (3.27a). 
For n — 3, the relevant family consists of three elements. 



the equation has the form 



C 12 
2 C 
6 3 



C 



(B22) 



with the solutions given in Eq. ( |6.27b ). 

For n = 4, the relevant family consists of five elements. 



the equation has the form 
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(B23) 
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with the solutions given in Eq. ( |6.27c ). 

For n = 5, the relevant family consists of seven elements, 



the equation has the form 
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with the solutions given in Eq. ( |6.27d ). 

For n = 6, the relevant family consists of eleven elements, 



the equation has the form 



= 



(B24) 
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with the solutions given in Eq. ( |6.27^ ). 
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TABLE I. Canonical dimensions of the fields and parameters in the model (|2.l| 
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FIG. 2. Leading scaling exponents for the sectors Z = 2, 8 and 12 (from above to below) in d = 2 (left) and d = 3 (right). 
Dashed lines denote solutions that exist as limits d — » 2 but disappear in two dimensions. 
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FIG. 3. Critical dimensions Ai-Ae (from below to above) of the operators (6.7) in the order 0(e) vs the space dimensionality 
d for 2 < d < 5 (left) and 4 < ci < 30 (right). The empty circles denote the operators which become trivial in d = 2 and 3. The 
dimensions tend to (Ai,2), e (A3-A5) and 2e (Ag) for d — > 00. 




in the scalar (solid 



curve) and vector (thick dots) models. Dashed lines denote monotonous branches of the critical dimensions in the vector case. 
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